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Abstract 

A quantum walk is a time-homogeneous quantum-mechanical process on a graph defined 
by analogy to classical random walk. The quantum walker is a particle that moves from a 
given vertex to adjacent vertices in quantum superposition. Here we consider a generaliza- 
tion of quantum walk to systems with more than one walker. A continuous-time multi-particle 
quantum walk is generated by a time-independent Hamiltonian with a term corresponding to a 
single-particle quantum walk for each particle, along with an interaction term. Multi-particle 
quantum walk includes a broad class of interacting many-body systems such as the Bose- 
Hubbard model and systems of fermions or distinguishable particles with nearest-neighbor 
interactions. We show that multi-paiticle quantum walk is capable of universal quantum com- 
putation. Since it is also possible to efficiently simulate a multi-particle quantum walk of the 
type we consider using a universal quantum computer, this model exactly captures the power 
of quantum computation. In principle our construction could be used as an architecture for 
building a scalable quantum computer with no need for time-dependent control. 



Introduction 



Quantum walk is a versatile and intuitive framework for developing quantum algorithms. Appli- 
cations of continuous- (i) and discrete-time (2, 3) models of quantum walk include an example of 
exponential speedup over classical computation (4) and optimal algorithms for element distinct- 
ness (5) and formula evaluation (6). 

Quantum walk can also be viewed as a model of computation. From this perspective it is 
natural to ask which quantum computations can be performed efficiently by quantum walk. This 
question was answered in reference (7) where it was shown that continuous -time quantum walk on 
sparse unweighted graphs is equivalent in power to the quantum circuit model. 

While this universality construction shows that quantum walk is a powerful computational 
model, it does not directly provide an architecture for a quantum computer. Since each vertex of 
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the underlying graph corresponds to a basis state in the Hilbert space, the graph associated with a 
quantum computation on n qubits is exponentially large as a function of n. This means that this 
quantum walk cannot be efficiently implemented using an architecture where each vertex of the 
graph occupies a different spatial location. Although the quantum walk on any sufficiently sparse 
graph can be efficiently simulated by a universal quantum computer (4, 8-10), it may be no easier 
to implement this quantum walk than to perform a general quantum computation. 

Nevertheless, many experimental implementations of quantum walk have been carried out. 
Typically, such experiments rely on non-scalable encodings of the quantum walk into a physical 
system such that locality of the graph translates into locality of the Hamiltonian, but with sub- 
stantial overhead that precludes the possibility of dramatic computational speedup (see for exam- 
ple (11-13)). While such experiments may serve as useful testbeds of presently available quantum 
information processors, they cannot directly realize efficient universal quantum computation. 

In this paper we consider a natural generalization of quantum walk. Whereas a standard quan- 
tum walk concems a single walker moving (in superposition) on a graph, we consider the general- 
ization to many interacting walkers. We show that such a multi-particle quantum walk is universal 
for quantum computation. Specifically, we show that any n-qubit circuit with g gates can be simu- 
lated by the dynamics of 0{n) particles that interact for a time poly(n, g) on an unweighted planar 
graph of maximum degree 4 with poly(n, g) vertices. Our construction is based on a discrete ver- 
sion of scattering theory and implements single-qubit gates similarly to the single-particle univer- 
sality construction (7). However, we use a different encoding of quantum data, and we implement 
two-qubit gates via interactions between pairs of particles. For indistinguishable particles, almost 
any interaction gives rise to universality. We present explicit universal constructions based on the 
Bose-Hubbard model, fermions with nearest-neighbor interactions, and distinguishable particles 
with nearest-neighbor interactions. 

To prove that our construction works, we develop several tools for analyzing multi-particle 
scattering on graphs. We prove error bounds for the propagation of one- and two-particle wave 
packets, and we prove a truncation lemma formalizing the idea that a particle moving at a fixed 
speed only "sees" part of the graph in a fixed time interval. Using these tools we prove that 
multi-particle quantum walk can efficiently simulate a standard quantum computer, establishing 
the computational power of multi-particle quantum walk. This result also provides limitations 
on the classical simulation of many-body interacting systems. For example, assuming quantum 
computers are more powerful than classical ones, our work implies that the dynamics of the Bose- 
Hubbard model on a sparse, planar graph cannot be efficiently simulated on a classical computer. 

Prior to this work it was known that systems of interacting particles on a lattice can be used 
to perform universal computation by changing the Hamiltonian as a function of time (14, 15). 
Because of the time dependence, such a system is not a quantum walk. In contrast, in our scheme 
the Hamiltonian is time-independent and the computation is encoded entirely in the graph on which 
the particles interact. 

Multi-particle quantum walk has previously been considered as an algorithmic tool for solving 
graph isomorphism (see for example (16)). However, this technique has known limitations (see for 
example (17)). Other previous work on multi-particle quantum walk has focused on two-particle 
quantum walk (18-26) and multi-particle quantum walk without interactions (18-22, 24, 27). Here 
we consider multi-particle quantum walks with interactions, which seem to be required in order 
to achieve efficient computational universality. Although non-interacting bosonic quantum walks 
may be difficult to simulate classically {28), such systems are probably not capable of perform- 
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ing universal quantum computation. For fermionic systems the situation is even clearer: non- 
interacting fermions can be efficiently simulated on a classical computer (29). 

We hope that our work motivates further experimental investigations of multi-particle quan- 
tum walk. Because our graphs are exponentially smaller (as a function of n) than those used in 
the single-particle construction, the multi-particle quantum walks we describe can be efficiently 
implemented using an architecture where vertices of the graph are represented by devices at differ- 
ent spatial locations (although we have not addressed issues of fault tolerance). The two-particle 
bosonic quantum walk experiments of references (20-22, 24) can be viewed as a first step toward 
implementing our construction. However, these experiments only involve non-interacting particles; 
as we discussed earlier, a nontrivial interaction appears necessary for universality. Conversely, our 
work shows that almost any interaction can be used to perform efficient universal computation. A 
Bose-Hubbard model of the type we consider could naturally be realized in a variety of experi- 
mental systems, including traditional nonlinear optics {30), neutral atoms in optical lattices {31) 
(with a lattice pattem implemented using a quantum gas microscope {32)), or photons in arrays of 
superconducting qubits {33). 



Multi-particle quantum walk 

In a multi-particle quantum walk, the particles interact in a local manner on a given simple graph 
G with vertex set V{G) and edge set E{G). We consider quantum walks with distinguishable or 
indistinguishable particles, where in the latter case the particles can be either bosons or fermions. 
The Hilbert space for m distinguishable particles on G is spanned by the basis 

G V{G)} (1) 

where iw is the location of the wih. particle. A continuous-time multi-particle quantum walk of m 
distinguishable particles on G is generated by a time-independent Hamiltonian 

t«=i (i,i)ei?(G) ^ ^ i,iey(G) 

where the subscript w indicates that the operator acts on the location register for the wth particle 
(tensored with the identity on all other particles). Here Uij{hi, hj) is a function of the number 
operators fii and hj that count the numbers of particles located at vertices i and j, respectively 
(explicitly, fn = X;™=i K)(^U)- 

The first term of (2) moves particles between adjacent sites, while the second term is an inter- 
action between particles. We only consider interaction terms acting between two or more particles, 
so we assume that Uij is zero whenever one of its arguments evaluates to zero and that Un is zero 
if there is only one particle at vertex i. We also assume that the interaction Uij has a constant range 
G E {0, 1, 2, . . .} (i.e., Uij = whenever the shortest path between vertices i and j has more than 
G edges). In our universality construction we consider graphs G that include long paths. Each 
vertex in one of these paths has degree 2 in the full graph G, i.e., the path is only connected at 
either end to other vertices in G. We require that for vertices i and j on such a path, Uij depends 
only on the distance between i and j (a kind of translation invariance). Finally, we assume that the 
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Figure 1 : The graph G obtained by attaching semi-infinite paths to a graph G. We label the 
vertices on the semi-infinite paths as {x,j), with j E {1, . . . , A^} indexing the path and x G N = 
{1, 2, 3, . . .} denoting the distance from G. The vertices labeled (1, 1), (1,2), . . . , (1, A^) are 
vertices of the original graph G. 

norm of each term Uij is upper bounded by a polynomial in m. 

Note that for any U satisfying our assumptions, the Hamiltonian (2) for a single particle reduces 

to 

(i,j)6E(G) 

the Hamiltonian for a standard continuous-time quantum walk (namely, the adjacency matrix of 
G). 

States representing m indistinguishable particles can be represented in the basis (1) as states 
that are either symmetric (if the particles are bosons) or antisymmetric (if the particles are fermions) 
under the interchange of any two particles. Since the Hamiltonian (2) is symmetric, it preserves 
both the symmetric and antisymmetric subspaces. Restricted to the appropriate subspace, the 
Hamiltonian (2) generates a quantum walk of m bosons or m fermions on G. 

This framework includes several well-known interacting many-body systems defined on graphs. 
For example, it includes the Bose-Hubbard model, where the interaction term is Uij{hi,hj) = 
{U /2)5ijhi{ni — 1). It also includes systems with nearest-neighbor interactions, such as the model 
with interaction term nj) = U5{^i^j)^E(G)nihj. 

Single-particle and two-particle scattering 

In our scheme for universal quantum computation we design a multi-particle quantum walk where 
the dynamics can be understood by considering scattering events involving one or two particles. In 
this section we show how single-particle and two-particle quantum walks on certain graphs can be 
analyzed using a discrete version of scattering theory. 

Single-particle scattering Consider a single-particle quantum walk on an infinite graph G ob- 
tained by attaching a semi-infinite path to each of N chosen vertices of an arbitrary (A^ + m) -vertex 
graph G as depicted in Figure 1 . Here we discuss the single-particle scattering process where a 
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particle is initially prepared in a state that moves (under Schrodinger time evolution) toward the 
subgraph G along one of the semi-infinite paths. After scattering through the subgraph, the particle 
moves away from G in superposition along the semi-infinite paths. To understand this scattering 
process we discuss the scattering eigenstates of the Hamiltonian Hq\ 

Scattering states are states with definite incoming momentum. Given the graph G, for each 
momentum k E (— tt, 0) and each path j E {1,2,..., A^}, there is a scattering state |scj(/c)) with 
amplitudes 

{x, q\scj{k)) = e-'^^Sgj + e'^^'Sgjik) (4) 

on the semi-infinite paths (with the labeling (x, q) of vertices on the paths as in Figure 1), where 
the N X N matrix S{k) appearing in the above equation is a unitary matrix called the S-matrix. 
In Appendix A we show how to obtain the S-matrix and the amplitudes of \scj(k)) within the 
graph G (34). The state \sc j{k)) is an eigenstate of H^^ with energy 2 cos A;. We also show in 
Appendix A that the states \sc j{k)) are delta-function normalized. 

A wave packet is a normalized state with most of its amplitude on scattering states with momen- 
tum close to some particular value. The scattering state \sc j{k)) gives us information about how a 
wave packet with momentum near k located on the semi-infinite path j scatters from the graph G. 
The wave packet initially moves toward the graph with speed | ^ | = |2 sin A;| . After scattering, the 
wave packet moves away from G along each of the semi-infinite paths (in superposition), and the 
amplitude associated with finding the wave packet on path q is Sqj{k). 

la Appendix E we discuss scattering of single-particle wave packets in more detail. In partic- 
ular, we show that this picture of the scattering process is valid when the finite extent of the wave 
packets is taken into account, and when the infinite paths are truncated to be long but finite. 

Two-particle scattering Now consider two-particle scattering on an infinite path. Translation 
symmetry makes this system easier to analyze than more general two-particle quantum walks (see 
for example references (35, 36)). 

We now discuss the scattering of two indistinguishable particles initially prepared in spatially 
separated wave packets moving toward each other along a path with momenta ki E (— tt, 0) and 
k2 E (0, tt) (we discuss distinguishable particles in Appendix C.2). Due to conservation of en- 
ergy and momentum, the state of this system after scattering still consists of one particle with 
momentum ki and another with momentum k2, but now moving apart. Since the particles are in- 
distinguishable, there is no distinction between transmission and reflection of the particles, so the 
effect of the interaction is to change the global phase of the wave function after scattering (relative 
to the case with no interaction). We show in the Appendices how to calculate this phase given the 
two momenta, the interaction U, and the particle type (bosons or fermions). In Appendix E we 
show that this picture of the two-particle scattering process holds on a long (but finite) path and 
with finite wave packets. 

Computation by multi-particle quantum walk 

In this section we describe our scheme for performing quantum computation using a multi-particle 
quantum walk of indistinguishable particles (we present a refinement of our scheme that uses 
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distinguishable particles in Appendix C.2). We encode n logical qubits into the locations of n 
particles on a graph of size poly(n). In addition to these n particles, we use another particle to 
encode an ancilla qubit that facilitates two-qubit gates. We call the original n qubits computational 
qubits, and we call the ancilla qubit the mediator qubit. Time evolution of a simple initial state 
with the Hamiltonian corresponding to a suitably chosen graph G implements a quantum circuit 
on this encoded data. 

The quantum circuit to be simulated (on the n computational qubits) is given as a product of 
one- and two-qubit gates from a universal gate set consisting of single-qubit gates along with the 
controlled phase gate 
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Note that a controlled phase between computational qubits i and j can be expanded as the following 
set of gates also acting on the mediator qubit m (initialized to |0)): 

CP,,|a„6„0^) = CNOT,^CP,^CNOT,^|a„6„0^) 

= HmCP^j^HmCPjmHinCP^jjjHinlaj, bj, Om). 

Writing the controlled phase gate in this way, we can transform the given quantum circuit into 
one with only single-qubit gates on computational qubits, Hadamard gates on the mediator qubit, 
and controlled phase gates between the mediator qubit and arbitrary computational qubits. In our 
construction we use this gate set acting on n + 1 qubits to simulate a quantum computation on n 
qubits. 

Each of the n + 1 qubits is represented in a dual-rail encoding using two paths that run through 
the graph, as shown in Figure 2. The encoded state |0) has a particle moving along the top path 
whereas the encoded state |1) has a particle moving along the bottom path. The particle moves as 
a wave packet with momentum near k. For the n computational qubits we choose k = — 7r/4 and 
for the mediator qubit we choose k = — 7r/2 (for concreteness). 

To implement single-qubit unitaries on each of the encoded qubits, we design the graph so 
that the particles scatter through a series of small subgraphs while remaining far apart. When the 
particles are all far from each other, the interaction term in the Hamiltonian is negligible and the 
n+1 wave packets propagate independently through the graph (here we neglect the interaction only 
for ease of exposition; our detailed analysis in Appendix D handles the full interacting system). In 
this case the multi-particle quantum walk can be viewed as a single-particle quantum walk for each 
of the particles. 

To implement the controlled phase gate between the mediator and a computational qubit, we 
design a subgraph that routes two particles toward each other and causes them to interact on a long 
path for a short time. The two particles then scatter away from one another and the system returns 
to a state where the particles are all far apart. 

Single-qubit gates To apply a one-qubit unitary U to a computational qubit, we insert an associ- 
ated graph G into the paths representing the qubit as follows. We attach two long "input" paths and 
two long "output" paths to four suitably chosen vertices of G so that the S-matrix at momentum 
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(a) Encoded |0). (b) Encoded |1) 



Figure 2: A qubit is encoded using single-particle wave packets. For an n-qubit computation we 
use n particles with momentum k = —n/A and one with k = —n/2. 




Figure 3: One-qubit gates at A; = — 7r/4. (a) Phase gate, (b) Basis-changing gate. 



— TT /4 has the form 

.^(«-). 

Each block of this matrix has size 2 x 2. A particle incident on the input paths with momentum 
— 7r/4 transmits perfectly to the output paths, with amplitudes determined by the unitary U. The 
scattering process implements the unitary U on the encoded qubit. 

Graphs G that implement a phase gate and a basis-changing gate at momentum — 7r/4 are 
shown in Figure 3 (7). The input and output paths are attached to the vertices denoted by open 
circles. The S-matrix at momentum — 7r/4 for each of these graphs is a 4 x 4 matrix of the form 
(5), with the lower left 2x2 submatrix given by 

t^phase Q t^basis ^ \-i I J' 

These two gates allow us to approximate arbitrary single-qubit unitaries on each of the n com- 
putational qubits. To obtain a graph implementing two of these unitaries in series, simply con- 
catenate the two graphs by attaching the output vertices of the first graph to the input vertices of 
the second graph. At momentum — 7r/4, the S-matrix of the resulting graph is the product of the 




Figure 4: Graph implementing a Hadamard gate at = — 7r/2. 
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S-matrices of the two graphs being concatenated (this is true whenever the two graphs have perfect 
transmission from input vertices to output vertices (7)). 

The Hadamard gate is the only nontrivial single-qubit gate we apply to the mediator qubit. To 
do so we use the graph depicted Figure 4 (37). This graph has S-matrix at momentum k = —n/2 
of the form (5), with lower left submatrix 

which is the Hadamard gate up to an irrelevant global phase. 

Two-qubit gates To implement the controlled phase gate between the mediator qubit and a com- 
putational qubit we use some facts about two-particle scattering on a long path. Recall that two 
indistinguishable particles of momentum ki and k2 initially traveling toward each other will, af- 
ter scattering, continue to travel as if no interaction occurred, except that the phase of the wave 
function is modified by the interaction. In general this phase depends on ki and k2 (as well as the 
interaction U and the particle statistics). For us, ki = —it/2 and k2 = tt /A (moving in opposite 
directions). We write e*^ for the phase acquired at these momenta. 

In our scheme we design a subgraph that routes a computational particle and a mediator particle 
toward each other along a long path only when the two associated qubits are in state |11). This 
allows us to implement the two-qubit gate 
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For some models CO = CP. We show in Appendix B that this holds in the Bose-Hubbard model 
(where the interaction term is V(ij{ni,nj) = {U /2)5ijhi{hi — 1)) when the interaction strength 
is chosen to be f/ = 2 + \/2, since in this case e*^ = —i. For nearest-neighbor interactions 
with fermions, with Uij{hi,hj) = U5(^ij-j^E{G)^i'ri'j, the choice U = —2 — \f2 gives e*^ = z, 
so CP = (C6')^. While tuning the interaction strength makes the CP gate easier to implement, 
almost any interaction between indistinguishable particles allows for universal computation. We 
can approximate the required CP gate by repeating the C6' gate a times, where e*'*^ ~ —% (which is 
possible for most values of Q, assuming Q is known (35)). 

Our strategy requires routing the particles onto a long path. This is done via a subgraph we call 
the momentum switch, as depicted in Figure 5(a). The S-matrices for this graph at momenta — 7r/4 
and — TT /2 are 

/ e-*^/^\ /I \ 

^switch (-7r/4) = 0-10 ^switch (-7r/2) = -1 . (6) 

\e-^'^/^ 0/ \0 -1 / 

The momentum switch has perfect transmission between vertices 1 and 3 at momentum — 7r/4 and 
perfect transmission between vertices 2 and 3 at momentum — 7r/2. In other words, in the schematic 
shown in Figure 5(a), the path a particle follows through the switch depends on its momentum. A 
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Figure 5: (a) Momentum switch, (b) C6 gate. 



particle with momentum —n/2 follows the double line, while a particle with momentum — 7r/4 
follows the single line. 

The graph used to implement the CO gate has the form shown in Figure 5(b). We specify the 
number of vertices on each of the paths in Appendix D. To see why this graph implements a C^^ 
gate, consider the movement of two particles as they pass through the graph. If either particle 
begins in the state |Oin), then it travels along a path to the output without interacting with the 
second particle. When the computational particle (qubit c in the figure) begins in the state |lin)'^, 
it is routed downward as it passes through the top momentum switch (following the single line). 
It travels down the vertical path and then is routed to the right (along the single line) as it passes 
through the bottom switch. Similarly, when the mediator particle begins in the state |lin)™, it is 
routed upward (along the double line) through the vertical path at the bottom switch and then to 
the right (along the double line) at the top switch. If both particles begin in the state |lin), then 
they interact on the vertical path. In this case, as the two particles move past each other, the wave 
function acquires a phase e*^ arising from this interaction. 

Note that timing is important: the wave packets of the two particles must be on the vertical 
path at the same time. We achieve this by choosing the number of vertices on each of the segments 
in the graph appropriately, taking into account the different propagation speeds of the two wave 
packets (see Appendix D for details). 

Discussion Our scheme for simulating a quantum circuit is summarized as follows. We encode 
n computational qubits as single-particle wave packets with momentum — 7r/4 traveling along 
two paths, along with a single mediator qubit similarly encoded but with momentum —n/2. We 
perform single-qubit gates as shown in Figure 3 and Figure 4, and we implement two-qubit gates 
using the mediator qubit and the graph shown in Figure 5. The subgraphs representing circuit 
elements are connected by paths. To illustrate how these ingredients are put together, a graph 
corresponding to a simple 2-qubit computation is depicted in Figure 6. 

In Appendix D we present all the details of our scheme and we prove that it performs the 
desired quantum computation up to an error term that can be made arbitrarily small. To prove this 
error bound we analyze the full {n + 1) -particle interacting many -body system. Our analysis goes 
beyond the scattering theory discussion presented in the previous section; we take into account the 
fact that the wave packets are finite (and therefore have support on a range of momenta), as well as 
the fact that the graphs involved in our scheme are finite. Specifically, we prove that by choosing 
the size of the wave packets, the number of vertices in the graph, and the total evolution time to be 
polynomial functions of both n and g, the error in simulating an n-qubit, g-gate quantum circuit is 
bounded above by an arbitrarily small constant (Appendix E). The bounds we prove are sufficient 
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Figure 6: Schematic depiction of a graph simulating B2 CPi 2 Ti on two qubits, for the Bose- 
Hubbard model with interaction strength [7 = 2 + \f2. The dotted lines represent paths and the 
single-qubit unitary gates represent their corresponding subgraphs {B is the basis-changing gate, 
T is the phase gate, and H is the Hadamard gate). 

to establish universality with only polynomial overhead, but they are almost certainly not optimal. 
For example, for the Bose-Hubbard model and for models with nearest-neighbor interactions, we 
prove that the error can be made arbitrarily small by choosing the size of the wave packets to be 
0{rP'g^\ the total number of vertices in the graph to be 0{rt}'^g^\ and the total evolution time to 
beC(niV)- 

We also describe two refinements of the scheme presented here. In Appendix C. 1 we show how 
the scheme presented above can be adapted so that the graph is planar and has maximum degree 
four, making it more amenable to implementation. In Appendix C.2 we show how the scheme can 
be adapted to use distinguishable particles with nearest-neighbor interactions. 

We have shown how multi-particle quantum walk can be used to perform efficient universal 
quantum computation. Our results provide an impetus for further experimental investigation of 
multi-particle quantum walk. We also hope that some of the tools we have developed in this work 
can be used to design and analyze new quantum algorithms based on multi-particle quantum walk. 
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A Single-particle scattering states 



In this section we establish some basic facts about the scattering states of the single-particle Hamil- 



tonian Hq \ Write 



H 



A 5t 
B D 



for the adjacency matrix of G, where A is N x N , D is m x m, and Bismx N . The first N rows 
and columns of this matrix correspond to the vertices attached to the semi-infinite paths. 

The results of this section apply to a broader class of Hamiltonians Hq^ than those considered 
so far. In particular, the results are valid when the finite graph G is a weighted, directed graph 
with edges in opposite directions having complex conjugate weights, in which case the Hamilto- 
nian H is an arbitrary Hermitian matrix. In this situation the full Hamiltonian Hq'^ has a term 
corresponding to H and a term representing the (unweighted, undirected) adjacency matrix of the 
semi-infinite paths. 

The S-matrix and the amplitudes of |scj(A;)) within the graph G can be obtained using the 
following procedure {34). Let il>j{k) be the m-component vector containing the amplitudes of 
\scj{k)) at the m vertices of G not connected to the infinite paths, and let 

Let z = e^^ and define 

zA - 1 zB^ \ 
zB zD-z^ -l) ' 



7(z) = 

Then the S-matrix and the amplitudes \I/ are given by 



f Sik) \ 1 _i 



(7) 



An equivalent expression for the S-matrix is 

Sik) = -Q{z)-'Qiz-') (8) 

where 



1-z(a + B^-. ^ b] . 



Using this expression, one can verify that the graphs described in this paper implement the claimed 
unitaries. However, it is not clear how to construct a graph that implements a desired unitary. In 
some cases one can find a suitable graph by performing a numerical search over graphs with a 
small number of vertices. Graphs with known properties at a given momentum can sometimes be 
combined (for example, the graph in Figure 8 was obtained in this manner). 
Using the expression 

7(2;)"^ = . . f adj7(z), 
det7(2;) 
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where adj 7(2) is the adjugate matrix of ■y{z), we see from (7) that each matrix element of S{k) and 
each matrix element of \l/(/c) is a rational function of z. Since S{k) is unitary, its matrix elements 
(as functions of z) have no poles on the unit circle. The following lemma shows that the matrix 
elements of ^(A;) also have no poles when z is on the unit circle. 

Lemma 1. Given G, there exists a constant A G M such that \{v\scj{k)) \ < X for all k G [— vr, tt), 
j G {1, . . .,N}, andv G G. 

Proof. Note that 
where 



0^ 

projects onto the N vertices of G attached to semi-infinite paths. Hence 



Let {{ipc) : c G {1, . . . , Uc}} be eigenstates of H satisfying -P|^c) = 0, and let this set be an 
orthonormal basis for the span of all such states. Then 

since each |^c) is an eigenvector of 7(2;) and Plipc) = 0. Reference (34) shows (in Part 2 of the 
proof of Theorem 1) that 

det ^(^) ^ for = 1' 

where M{z) is the {N + m — nc) x (N + m — nc) matrix of 7(2;) in the subspace of states orthogonal 
to the span of {\ipc) : c G {1, . . . , Uc}}. Therefore 

l{v\scj{k)) = -(t;|7(2)"S ( 7 ) \j) 



z \ z 



{v 



has no poles on the unit circle, and the result follows. □ 
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We now establish the delta-function normalization of the scattering states. Let 

oo N 

ni = ^^\x,q){x,q\ 

x=l q=l 

oo N 

^2 = ^-^^\x,q){x,q\ 

x=2 q=l 

N 

n3 = 

q=l 

We show that, for k e (-tt, 0), p e (-tt, 0), and i,j e {1,..., N}, 

{sci{p)\scj{k)) = {sCi{p)\Ui + Ila - U^lsc j{k)) = 27r5ij5{k - p). 

First write 

oo N 
x=l q=l 



^ / N \/oo oo \ 

2 + E sup)Sqjik) Yl + E ^"'^'"'^^ 

\ q=l J \x=l x=l J 

^ / N \ / oo oo 

+ 2 - E ^W^^Ak) E - E 

V q=l / \a:=l x=l 



^-i{p-k)x 



oo oo 

i(j>+k)x 



+ 2('^;.(^^)+'^^.(^)) E^"^^"^'^'+E 

\x=l x=l 
^ / oo oo \ 

+ ^is*,{p) - s,j{k)) Yl - E 

\x=l x=l J 

We use the following identities for p, A; e (— tt, 0) : 

oo oo 

^ gi(p-fe)a: ^ ^ g-i(p-fe)x ^ 27r5(p - A;) - 1 

x=l x=l 
oo oo 

gi(p+fe)a; _|_ _ _-|^ 

a;=l a;=l 

^ ^ /r)-k\ 

J2 e'^P-'^^ - J2 = i cot ( ^ J 

x=l x=l ^ ^ 

^ e^(^'+^)^ - Y = ^ cot . 

x=l x=l ^ ^ 
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These identities hold when both sides are integrated against a smooth function of p and k. Substi- 
tuting, we get 

(scj (p) I Hi I scj (k)) = 27r% 5(p - /c) + % ( ^ cot ' ^ 



where we used unitarity of the ^'-matrix to simplify the first term. Now turning to 112 we have 

{sCi{p)\H]l2\sCj{k)) = 2cos(p)(sCi(p)|n2|sCj-(/c)) 

and 

N 



2/2 

2 



{sc,{p)\HIl2\sc,{k)) = (sc,(p)| ( 2cos(fc)n2|sc,(A;)) + J^i'^'^'^d + S,j{k)e''')\'2, q) 

^ q=l 

N . 

q=l / 



Q 

Using these two equations we get 

(2cos(p) - 2cos(A;))(sCi(p)|n2|sc,(A;)) = 6ij{e^'P-"' - e-^^'+'P) + S*,{p){e-^'P-"' - e-^^^-'P) 

N 
q=l 



Noting that 



AT 

(sc,(p)|n3|sc,(^)) = J2(^iqe'P + S;,ip)e-'n{Sjqe-'' + S,,{k)e% 

q=l 



17 



we have 

2co.(ri"-2cos(^ --n 

1 1 -ip-ik \ 

2cos(p) -2cos(A;) J 

(^2ip+ik g2ifc+jp \ 
2cos(p)-2cos(A;) ~ ''''^') 

/ -2ip+ik _ 2ik-ip \ 

+ E S,'.(p)Sb(*) + ^ (^) ) ■ (11) 
Adding equation ( 10) to equation (11) gives equation (9). 



B Two-particle scattering states 

Here we derive scattering states of the two-particle quantum walk on an infinite path. We write 
the Hamiltonian in the basis \x,y), where x denotes the location of the first particle and y denotes 
the location of the second particle, with the understanding that bosonic states are symmetrized and 
fermionic states are antisymmetrized. The Hamiltonian (2) can be written as 

= Hi'^ ® I, + I, ® + 5^ V(|x - yl) \x, y) {x, y\ (12) 

where V corresponds to the interaction term U and (with a slight abuse of notation) the subscript 
indicates which variable is acted on. Here 

H^^^ = ^ |x + + \x){x + 1| 

is the adjacency matrix of an infinite path. Our assumption that U has finite range C means that 
V(r) = Oforr> C. 

The scattering states we are interested in provide information about the dynamics of two parti- 
cles initially prepared in spatially separated wave packets moving toward each other along the path 
with momenta ki E (— vr, 0) and k2 E (0, vr). 

We derive scattering eigenstates of this Hamiltonian by transforming to the new variables s = 
X + y and r = x — y and exploiting translation symmetry. Here the allowed values (s, r) range 
over the pairs of integers where either both are even or both are odd. Writing states in this basis as 
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s;r), the Hamiltonian takes the form 



® + I, ® ^ V(|r|) |r)(r 



(13) 



For each pi e (— tt, tt) and p2 e (0, tt) there is a scattering eigenstate |sc(pi; P2)) of the form 



where the state \ip{pi;p2)) can be viewed as an effective single-particle scattering state of the 
Hamiltonian 



with eigenvalue 4 cos(pi /2) 005(^2) • For a given V, the state | V' (Pi ; P2) ) can be obtained explicitly 
by solving a set of linear equations (see for example (34)). It has the form 



for P2 G (0, tt). Here the reflection and transmission coefficients R and T and the amplitudes of 
the scattering state for |r| < C (described by the function /) depend on both momenta as well as 
the interaction V. With R, T, and / chosen appropriately, the state |sc(pi;p2)) is an eigenstate of 
i?^^) with eigenvalue 4 cos (pi/2) cos(p2)- 

Since V(|r|) is an even function of r, we can also define scattering states for p2 e (— tt, 0) by 



These other states are obtained by swapping x and y, corresponding to interchanging the two 
particles. 



(s;r|sc(pi;p2)) = e-^^'^^/^(r|V^(pi;p2)), 




(14) 




(15) 



(s;r |sc(pi;p2)) = (s; -r|sc(pi; -^2))- 
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Figure 7: Scattering of two particles on an infinite path. 
The states {|sc(pi;]?2)) : Pi G (— tt, tt), p2 ^ (— tt, 0)U(0, tt)} are (delta-function) orthonormal: 
(sc(p'i;P2)|sc(pi;P2)) = {sc{p[; p'^)] I ^ |r)(r| ® J |sc(pi;p2)) 

\r, s even / 

+ (sc(pi;P2)| I ^ I |sc(Pi;P2)) 

\r, s odd / 

^ ^ e-^^'-^i)^/^ 5^(V^(p;;p'2)|r)(r|^(pi;p2)) 

s even r even 

+ J2^''^"~''^'^'J2WilP2)\r){rmPi,P2)) 

s odd r odd 

oo 

^27r6{pi-p[) (V'(pi;P2)|r)(r|V^(pi;p2)) 

r=— oo 

^ ^n''S{pi - p[)S{p2 - P2) 

where in the last step we used the fact that {iApiW^) P2)) = 2n5{p2 — p^)- 

To construct bosonic or fermionic scattering states, we symmetrize or antisymmetrize as fol- 
lows. Forpi e (— TT, tt) andp2 £ (0, vr), we define 

1 



|sc(pi;p2))± = --^(|sc(pi;p2)) ± |sc(pi; -^2))). 



Then 



' Q-ip2r _|_ gje±(pi,P2)gip2r- if r < — C 

(s;r|sc(pi;p2))± = -^e-^^'i^/2<j/(pi,P2,r)±/(pi,P2,-r) if|r| <C (16) 

^ I gi6'±(pi,P2)g-«P2r- _|_ gip2r- if r > C 



where 9±{pi,p2) is a real function defined through 

e''±^^'^^)=r(pi,P2)±i?(pi,P2). 



(17) 



Note that \T ± R\ = 1; this follows from the potential V(|r|) being even in r and from unitarity 
of the S-matrix. These eigenstates allow us to understand what happens when two particles with 
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momenta ki E (— vr, 0) and k2 G (0, tt) move toward each other. Here pi = —ki — k2 and 
P2 = {k2 — ki)/2. Recall (from the main text of the paper) that we defined e*^ to be the phase 
acquired by the two-particle wavefunction when ki = —tt/2 and ^2 = n/A (9 depends implicitly 
on the interaction V and the particle type), so 6* = 9±{n/4:, Svr/S) . 

For |r| > C the scattering state is a sum of two terms, one corresponding to the two particles 
moving toward each other and one corresponding to the two particles moving apart after scattering. 
The outgoing term has a phase of T ± i? relative to the incoming term (as depicted in Figure 7). 
This phase arises from the interaction between the two particles. 

For example, consider the Bose-Hubbard model, where V(|r|) = U5r,o- Here C = and 
T = 1 + R.ln this case the scattering state |sc(pi; ^2))+ is 

-1 

{x,y\sc{pi;p2))+ = -^e-*P^(^) (e^P^I-^/l + e^^+^^i-P^^e-'^'^l"-^!) . 

The first term describes the two particles moving toward each other and the second term describes 
them moving away from each other. To solve for the applied phase e^^+'^^^'P^^ we look at the 
eigenvalue equation for P2)) at r = 0. This gives 

U 



U — Ai cos(pi/2) sin(p2) 
So for the Bose-Hubbard model, 

i9+{pi,p2) _ r^r^ ^ \ , TD^^ ^ \ _ ^ + 4«cos(pi/2) sin(p2) _ 2 (sin(fc2) - sin(A;i)) - iU 



T{pi,P2) + R{Pl,P2) 



U — 4d cos(pi/2) sin(p2) 2 (sin(A;2) — sin(A;i)) + iU 



For example, if f/ = 2 + a/2 then two particles with momenta ki = —n/2 and k2 = 7r/4 acquire a 
phase of e^*'^/^ = —i after scattering. 

For a multi-particle quantum walk with nearest-neighbor interactions, V(|r|) = US\r\,i and 
C = 1. In this case the eigenvalue equations for P2)) at r = — 1, r = 1, and r = are 

4cos cos(p2)(e'^^ + R{pi,p2)e-'P') = U{e'P' + R{p^,p2)e-'P') 

+ 2cos (I) {e''^- + R{p,,p2)e-''P-' + f{puP2,0)) 
4cos (1^) cos{p2)T{p,,p2)e-'^' = UT{p^,p2)e-'''' 

+ 2 cos (1^) (/(pi,P2,0)+T(pi,p2)e-2^^^) 
2cos(p2)/(pi,P2,0) = T{puP2)e-'^' + e'"' + R{pi,P2)e-'^\ 

respectively. 

Solving these equations for R, T, and f{pi,P2, 0), we can construct the corresponding scat- 
tering states for bosons, fermions, or distinguishable particles (for more on the last case, see 
Appendix C.2). Unlike the case of the Bose-Hubbard model, we may not have 1 + R = T. 
For example, when U = —2 — \f2, pi = 7r/4, and p2 = Stt/S, we get i? = and T = i (see 
Appendix C.2). 
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in 



out 



Figure 8: A planar graph that implements a Hadamard gate at momentum —-n/2. 

C Refinements of the universafity construction 

In this section we present two refinements of our scheme. We first show how our scheme can be 
modified to use planar graphs of maximum degree four. We then give a universality construction 
using distinguishable particles with nearest-neighbor interactions. 

C.l Making the graph planar 

The example in Figure 6 shows that the graphs in the scheme described by the main text of our 
paper may not be planar: the mediator qubit can interact with any of the computational qubits, so 
the vertical paths for the CO gate cross other paths in the graph. Furthermore, both the graph used 
to implement the Hadamard gate on the mediator qubit (Figure 4) and the graph used to implement 
the gate (Figure 5(b)) can lead to a nonplanar overall graph when input and output paths are 
attached in the prescribed manner. In this section we describe how to modify the scheme so that 
the resulting graph is planar and has maximum degree 4. 

The first simple modification is to replace the graph from Figure 4 with a planar graph (with 
input and output vertices on the same face in the correct relative positions) that also implements 
a Hadamard gate on the mediator qubit. The graph in Figure 8 does the trick: its S-matrix at 
momentum — 7r/2 has the form (5) with lower left submatrix 



The smaller maximum degree (4 instead of 5) and planarity of this graph come at the expense of 
increasing the number of vertices (as compared to the graph in Figure 4). 

As a second modification, we introduce additional mediator qubits. Throughout the graph, 
we arrange the input and output paths for the computational qubits vertically from 1, . . . , n with 
the path corresponding to logical always above the path corresponding to logical 1. For each 
i E {1,...,?7, — l}we place a mediator qubit labeled m(i) between computational qubits i and 
i + 1. We only perform two-qubit gates between adjacent qubits throughout the computation (i.e., 
mediator qubit m(z) only interacts with logical qubits i and i + 1). 

To implement two-qubit gates in a planar manner, we use the graph shown in Figure 9. This 
graph is obtained by concatenating two CO graphs and uncrossing paths to make the drawing 
planar. We only use this gate between adjacent encoded qubits, one of which is a mediator qubit 
and one of which is a computational qubit. Note that this graph involves two adjacent paths (path 
1 of the top encoded qubit and path of the bottom encoded qubit) as opposed to the two 1 
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li+l.in °- 



<^ 



^ lj+l,out 



Figure 9: The planar entangling gate between adjacent encoded qubits m(z) and i + 1. This 
graph implements the unitary Xj+i(C6')^j^^^ ^^-,^Xj+i. A similar graph implements the unitary 
Xm(i)(C6')^jj^(-)Xm(j) between adjacent encoded qubits i and m{i). 

paths in the C9 gate in Figure 5. The resulting logical gate is (C^^)^ conjugated by an X gate on 
the bottom qubit. More explicitly, if we interact computational qubit i and mediator ra{i), we 
implement Xin(j)(C^)^j^(.-)Xni(i), whereas if we interact mediator m(i) and computational qubit 
z + 1, we implement Xi+i(C6')^^.^ .^^Xj+i. Applying this gate a times, where e*^'* ^ ±i, we can 
approximate the gates Xm(i)CZi,m(i)Xm(i) and Xi+iCZm(i),i+iXi+i. 

Using these gates, we can perform a controUed-Z gate between computational qubits i and i + 1 
by the following sequence: 

CZj^j+i|aj, Om(i)) = CNOTj+i^in(i)CZj^ni(j)CNOTj+i^jn(j)|aj, 0^(1)) 

= Xjn(i)Xi+iCNOTj+i^ni(i)CZj^in(i)CNOTj_|_i^m(j)Xin(j)Xj+i|aj, Om(j)) 
= Hni{j) (Xj+iCZni(j),i+lXj+i) Hjn{j) (Xi„(j)CZj^i„(j)Xjn{j)) Hin(j) 
(Xj+iCZin(i),i+lXj+i) Hm(j)|aj, Oin(i))- 

To implement a CZ gate between arbitrary encoded qubits, we use these CZ^ i+i and one-qubit 
gates to implement a SWAPj j+i gate, facilitating movement of encoded qubits. To implement a 
CZj j gate, we use SWAP gates to move the information encoded in qubit i to qubit j — 1 or j + 1, 
perform the required CZ gate, and finally use SWAP gates to return qubit i to its original position. 

Note that every graph implementing a gate preserves the relative position of each path for both 
the input and the output (e.g., qubit 1 remains above qubit m(l) and each path remains above each 
1 path), so concatenation of the graphs preserves planarity. As each individual graph is planar, and 
as concatenating two graphs preserves planarity, the overall graph simulating the quantum circuit 
is planar. 

C.2 Distinguishable particles 

So far we have focused on the case of indistinguishable particles. However, we can also perform 
universal quantum computation with distinguishable particles, provided the interaction has an ap- 
propriate form. 

For distinguishable particles we use the same encoding of qubits as before (computational 
qubits have momentum —tc/A and mediator(s) have momentum — 7r/2), except that now each 
qubit is associated with a specific particle (e.g., computational qubit 1 is associated with particle 



23 



1). Since the graphs implementing single-qubit gates only make use of single-particle scattering, 
which is unaffected by particle statistics, we can use the same graphs for distinguishable particles. 
As such, we need only examine the implementation of the C^^ gate to see how our construction 
must be modified. In this section we show that with a simple nearest-neighbor interaction we 
can make our scheme work for distinguishable particles by carefully choosing the strength of the 
interaction term in the Hamiltonian, but with no other modifications. 

When two indistinguishable particles of momenta ki and ^2 scatter on an infinite path, there 
is no distinction between the final state where the particles reflect off of each other (exchanging 
momenta) and where the particles transmit through one another. Thus, after scattering, the global 
phase of the wave function is multiplied by a factor T ± R, the sum of the amplitude to transmit 
and the amplitude to reflect (or the difference if the particles are fermions). For any interaction 
potential, \T ± R\ = 1, and in most cases the applied phase is nontrivial and can be used for 
universal computation within our scheme. 

In contrast, when two distinguishable particles of momenta ki and k2 scatter on an infinite 
path, there are two distinct outgoing states: one corresponding to the case where the two particles 
reflect and one where the particles transmit. We circumvent this potential problem by choosing the 
interaction strength so that the transmission probability for two-particle scattering at momenta 7r/4 
and — 7r/2 is 1 (forcing R = 0), yet T 7^ 1 (so that T is a nontrivial phase). With such a choice, 
the graph implementing the C9 gates preserves our encoding of qubits. In other words, if encoded 
qubit 1 is associated with particle 1 before applying the gate, then it is still associated with particle 
1 after applying the gate (and similarly for the second qubit involved in the gate). We can then use 
the same graph as before to implement the controlled phase gate between encoded qubits. 

Consider the nearest-neighbor Hamiltonian with Wjj(nj, Wj) = U5ij(zE{G)ninj. For two parti- 
cles on an infinite path, this is (12) with V(|r|) = US\r\,i- The reflection coefficient R(pi,p2) for 
Pi = 7r/4 andp2 = 37r/8 is 



Our goal is to choose U so that R = and T is a nontrivial phase. The values of U that set 
R = are. U = Oorf/=— 2 — \/2. The solution U = corresponds to no interaction and the 
trivial phase T = 1 which is not sufficient for universal computation within our scheme. Choosing 
U = —2 — a/2 sets T = i which allows us to perform a CP gate. 

We expect that other types of multi-particle quantum walk with distinguishable particles can 
also be used for universal computation. However, unlike in the case of indistinguishable particles, 
the interaction term may have to be tuned to satisfy the conditions R = and T 7^ 1 as was the 
case here. For some interactions, it may not be possible to satisfy these two requirements. For 
example, for a model where interactions only occur when particles occupy the same site (such as 
in the Bose-Hubbard model) we have 1 + R = T for all momenta, so the transmission amplitude 
is trivial whenever R = 0. 

Note that it may be possible to implement an entangling two-qubit gate in other ways. For 
example, some interactions may allow two-particle scattering with T = and R = i, in which 
case the graph shown in Figure 9 preserves the encoding of qubits and implements such a gate. 
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D Detailed description of the scheme 



In this section we fill in all of the details of the scheme outlined in the main text of this paper. 
We specify the initial state, the graph used to perform an n-qubit quantum computation, and the 
evolution time. These are all specified as a function of a single parameter L E N. We show that 
by taking L = poly(r;,, g) we can achieve an arbitrarily small error in our simulation of a given 
g-gate quantum circuit. The resulting graph and evolution time are both polynomially large in n 
and g. For simplicity, we discuss the case of indistinguishable particles and we concentrate on the 
nonplanar scheme using one mediator qubit, but similar results clearly apply to the planar scheme 
and for the case of distinguishable particles as discussed in Appendix C.l and Appendix C.2. 

D.l Building the graph 

The graph corresponding to a given circuit is built by piecing together subgraphs (which we call 
blocks) that implement gates. These blocks are of the forms shown in Figure 10 (which we call a 
block of type I) and Figure 1 1 (which we call a block of type II). 

A block of type I applies single-qubit gates f/i, . . . , f/„ to the computational qubits and a 
single-qubit gate V^^ to the mediator, all acting in parallel. The unitaries f/i, . . . , f/„ may be 
phase gates, basis-changing gates, or identity gates, up to a global phase. The unitary V^a is ei- 
ther the identity gate or the Hadamard gate, up to a global phase. The circles in Figure 10 are 
replaced by corresponding subgraphs implementing the single-qubit gates. To simplify our analy- 
sis in Appendix E.5, we replace the circles with subgraphs with the property that the shortest path 
between the two output vertices (or input vertices) of the subgraph is greater than C, the interac- 
tion range of the Hamiltonian. To implement the identity gate in this manner we use a graph that 
connects each input vertex to the corresponding output vertex by an edge. Similarly, this condition 
is automatically satisfied by the subgraph given in Figure 3(a) for the phase gate (for any C), and is 
satisfied by the single-qubit subgraphs for the basis-changing gate (Figure 3(b)) and the Hadamard 
gate on the mediator (Figure 4) when C E {0, 1, 2, 3}. If C > 4 this condition can be achieved by 
simply concatenating the graph implementing the desired unitary with two paths of length C + 1 
on the input as well as the output. This has the same effect as simply adding extra vertices on 
the input and output paths of the block, but in Figure 10 we include these extra vertices inside the 
circled subgraph. 

For a block of type I, the number of vertices on each input path for the mediator qubit is equal 
to the number of vertices on each of its output paths and is chosen to be 



where M(— 7r/2) = L. Similarly, the number of vertices on the input and output paths is the same 
for each computational qubit and is chosen to be 



2M(-7r/2) + L 



2M(-7r/4) + L 



where 
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2M(-f)+L 2M(-f) + L 




Computational 
Qubits 



2M(-f) + L 2M(-f) + L 




Mediator 
Qubit 



Figure 10: A block of type I. Here the single-qubit gates C/i, . . . , C/„ on the encoded computational 
qubits are either the phase gate, the basis-changing gate, or the identity (up to a global phase). The 
single-qubit gate on the mediator is either the identity or the Hadamard gate (again, up to a 
global phase). 
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1 1 1 






L + 2M{k) 



M{k) L M{k) 



Figure 12: A single-qubit gate U acts on an encoded qubit. The wave packet starts on the paths on 
the left-hand side of the figure, a distance M(k) from the ends of the paths. After time tj = 3L/2 
the logical gate has been applied and the wave packet has traveled a distance 2M{k) + L (up to 
error terms that are bounded as 0{L-^/^)). 

These lengths are chosen to be different to compensate for the fact that the computational and 
mediator particles move at different speeds. They are designed so that a wave packet of length L 
and momentum k on an input path initially a distance M{k) from the subgraph implementing the 
unitary will be found a distance M{k) from the graph on the output paths after time ti = 3L/2. 
This is illustrated in Figure 12. 

Blocks of type II are used to perform gates between a computational qubit and the mediator. 
The two-qubit gate involved in such a block, including the lengths of all the relevant paths, is shown 
in detail in Figure 15. To implement an identity operation on the other computational qubits, we 
simply attach the input to the output by a path with 4M(— 7r/4) + 6L + 4 vertices. A wave packet 
of momentum — 7r/4 and length L starting on an input path a distance M(— 7r/4) from the leftmost 
vertex will be found a distance approximately M(— 7r/4) from the rightmost vertex on an output 
path after a time tu = (5L + 2M(-7r/4))/v^. 

To compose blocks B and B', where B' is a block of type I, replace each of the vertices of the 
input paths for block B' with the rightmost 2M{k) + L vertices of the corresponding output paths 
of block B. If B' is a block of type II, replace only the leftmost 2M [k) + L vertices of the input 
paths of B' with the rightmost 2M{k) + L vertices of the output paths of B. A simple example is 
shown in Figure 13. 

This method of composing blocks simplifies our analysis, which proceeds by computing the 
evolution of the particles inside each block and then using a "truncation lemma" to bound the er- 
rors arising from the connection to other blocks on either side. In the example shown in Figure 13, 
the initial state of each particle is a wave packet of length L prepared on the input paths of block B 
a distance M{k) from the single-qubit gate subgraphs. At time ti the single-qubit gates have been 
applied and the wave packets encoding the logical state have propagated a distance approximately 
M{k) on the output paths of the first block. This output state for the first block coincides (approx- 
imately, and up to an irrelevant global phase) with an input logical state for the second block, as it 
is a distance approximately M {k) from the second set of single-qubit gates. Finally, after time 2ti 
the logical state is on the output paths of the second block. 

More generally, for a circuit with gi blocks of type I and gn blocks of type n, the total evolution 
time for the computation is 

T = giti + giitii = Q{gL), (18) 

where fif = 5fi + 5fii. 
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Figure 13: An example of combining blocks to build a circuit. Here the two blocks are both of 
type I. 



D.2 Initial state, final measurement, and error bound 

We prepare an initial state of n + 1 spatially separated wave packets, where each wave packet is 
localized on an input path of the first block of the graph. To simplify notation, we assume that 
the first block is of type I, taking the unitaries Ui and to be the identity if required. Label the 
vertices of the 2n input paths in the first block corresponding to computational qubits as \x, q)l^, 
where j e {1, . . . ,n} indexes the qubit, q E {0, 1} labels the computational basis state, and 
X e {1, . . . , 2M(— 7r/4) + L} denotes the position along the path. Similarly, label the vertices of 
the input paths for the mediator qubit as \x, q)^ where x G {1, . . . , 2M(— 7r/2) + L}. Here the 
vertices with a: = 1 are the rightmost vertices on the input paths. We define the logical input states 
for each of the qubits as 



Ml 



Vl 



E 



x,0)l 



11- V 
111/ 



x=m{-^)+1 



E 

x=M(-f)+l 



and the logical input states for the mediator as 



Mi 



E 



=M( 



0+1 



Ml 



11- — 
J-in/ — 



E 



x=M[-^)+l 



The system is initialized at time t — Oin the computational basis state |00 . . . 0), encoded as 

|V^(0)) = Sym(|Oi„)^..|Oi„nOi„)°^). 

Here Sym is a linear operator that symmetrizes if the particles are bosons and antisymmetrizes if 
the particles are fermions. It is defined by 

Sym(|ai)|a2) ■ ■ ■ |a^+i)) = ^ = (±l)'§"('')|a^(i))|a^(2)) . . . |a^(n+i)) 

V(n + 1)! 
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where the ± is + for bosons and — for fermions. 

We evolve the initial state for a time T according to the Schrodinger equation with the Hamil- 
tonian Hq^^^ on a graph G built by composing gi blocks of type I and gn blocks of type II as 
described in the previous section. The total evolution time T is given by equation (18). Labeling 
the vertices on the output paths of the final block (again assumed to be type I) as |x, g)out> where 
now vertices with x = 1 are the leftmost vertices on the output paths, the logical output states for 
the computational qubits are defined as 



A/(-f )+L 
x=M(-f)+l 



3 

out 



II \J 
I J-out/ 



71 



E 



with similar definitions for the mediator qubit. Letting Uc be the logical n-qubit unitary that the 
graph is intended to implement, the desired output state is 

|0) = Sym[ {z'z\..z-\Uc\m...Q)\zl,)K.^z:^n^,,,y 
\2e{o,i}" 



where z"^ is the zth bit of z. A final measurement of the encoded quantum state in the computational 
basis is performed by measuring the locations of the particles at the end of the time evolution. 
We prove in the next section that 



1^(0)) -e' 



0(gn 



H. 



(n+l) 



(19) 



where 7 is an irrelevant overall phase. For the Bose-Hubbard model and for the models with 
nearest-neighbor interactions that we consider, = 0{n'^), so by taking L = 0{n^'^g'^) 

we can make the error arbitrarily small. Using this bound on L, the total number of vertices 
required in our construction is 0{n^^g^) and the total evolution time is 0(71^"^ g^). Although these 
bounds are sufficient to establish universality with only polynomial overhead, we expect that they 
can be improved significantly. 



E Analysis of wave packet scattering 

In this section we prove the error bound (19) stated above. The proof relies on the following 
technical results. 

In Appendix F we prove the following theorem about single-particle wave packet scattering 
on an infinite graph of the form shown in Figure 1. The proof is based on a calculation from 
reference {6). 

Theorem 1. Let G be an {N + m) -vertex graph. Let G be a graph obtained from G by attaching 
semi-infinite paths to N of its vertices, as shown in Figure 1, and let S be the corresponding S- 
matrix. Let Hq^ be the quantum walk Hamiltonian of equation (3). Let k G (— vr, 0), M, L G N, 
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j G {1, . . . , A^}, and 

. M+L 

Let Co Z^e a constant independent of L. Then, for allO < t < cqL, 



where 



Wit)) 



with 



oo N 



— 2it cos k 



Vl 



^ J] (5gje-*'=^i?(x - [2tsinA;J) + Sgj{k)e'''''R{-x - [2tsmk\)) \x,q) 

x=l q=l 



R{1) 



1 ifl e {M + 1,M + 2,. . .,M + L} 
otherwise. 



The approximation \a^{t)) to the time-evolved state has two terms, one corresponding to the 
incoming wave packet and one corresponding to a superposition of outgoing wave packets. Before 
scattering, the wave packet is supported entirely on path j and the second term is zero. Likewise, 
after scattering the first term is zero and the particle is outgoing along the semi-infinite paths. As 
we can see from the above expression, the square wave packets involved move with speed 2|sin A;| 
as expected (and to a good approximation, movement occurs at discrete times). 

In Appendix G we prove the following result about two wave packets moving past each other 
on an infinite path. While the same proof applies to wave packets with other momenta, for con- 
creteness we assume that the momenta are ki = —n/2 and /c2 = 7r/4 (so pi = —ki — k2 = it /A 
andp2 = (^2 - ki)/2 = Svr/S). 

Theorem 2. Let H^'^^ be a two-particle Hamiltonian of the form (12) with interaction range at 
most C, i.e., V(|r|) = for all \r\ > C. Let 0±{pi,p2) be given by equation (17). Define 9 = 
0±(7r/4, 37r/8). Let L e N, let M e {C + 1,C + 2, . . .}, and define 



\Xz,k) 



1 



ikx 



x) 



x=z—L 



IV^(O)) = ^ (|X-M,-f )|XAf+L+l,f ) ± |XM+L+l,f )|X-M,-f )) • 



Let Co be a constant independent of L. Then, for all < t < cqL, we have 



where 



\a 



(t)) = ^a^y{t)\x,y), 



(20) 
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a:j;y{t) = ±ay^{t)forx ^ y, and, for x < y, 

a.,y{t) = -^e-^'' [e-'^''/^e'^y'^F{x,y,t)±e''e''''''^e-'^y'^F{y,x,t 
v 2L 

where 



(21) 



F{u,v,t) 



1 ifu~2[t\e{-M -L,...,-M -l}andv + 2 ^ e {M + 1, . . . , M + L} 



otherwise. 



Note that for t > to, where to is the smallest time that satisfies 

V2_ 



[toj >2M + 2L + 1, 



we have F{x, y,t) = for x < y. Thus, for t > to. 



ie 



= (lX-Af+2LtJ,-f)|XM+L+l-2[t/v^J,f ) ± IXA./+L+l-2[t/v^J,f)IX-Af+2LtJ,-f) 

(22) 

which describes two wave packets of length L moving apart. In particular, note that the prefactor 
includes an overall phase of e*^ arising from the interaction between the particles. 

Note that in Theorems 1 and 2, the big-O notation includes the dependence on M. In particular, 
the error bounds hold even if M is chosen to depend on L. When applying these Theorems, we 
choose M to be proportional to L. 

These results about wave packet scattering on infinite graphs are not sufficient to prove error 
bounds for our scheme, as the graphs we construct are finite. However, a wave packet propagates 
with a finite speed, so that if we evolve for a finite amount of time then we expect a wave packet to 
only "see" a small part of the graph. The next ingredient in our analysis is a "truncation lemma" 
that we use to make this precise. The lemma, which we prove in Appendix H, is as follows: 

Lemma 2 (Truncation Lemma). Let H be a Hamiltonian acting on a Hilbert space T-L and let 
|<I>) E 7i be a normalized state. Let K be a subspace ofH, let P be the projector onto K, and let 
H = PHP be the Hamiltonian within this subspace. Suppose that, for some T > 0, W E {H, H}, 
Nq G N, and 5 > 0, we have, for all < t < T, 



iWt 



\^) = \l{t)) + \e{t))with |||e(t))||<5 



and 



1 - P)W'\-i{t)) = Ofor all r e {0,l,...,No- 1}. 



Then, for all < t < T, 



'iHt „-iHt 



1$) 



< I ^^t^ + 2) (5 + 2-^«(l + 5)). 
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E.l Truncating the semi-infinite paths 



Lemma 2 lets us prove analogs of Theorem 1 and Theorem 2 when the graphs involved have been 
truncated to have finitely many vertices. 

For example, consider the case where H = Hq^ is the Hamiltonian (3) for a single particle on 
a graph of the form shown in Figure 1 . Let G{K)hQ the finite graph obtained from G by truncating 
each of the paths to have a total length K = il{L) (so that the endpoints of the paths are labeled 
{K, j) for j G {1, . . . , A^}), and choose H = H^l^j^y Let the subspace /C be spanned by basis states 
corresponding to vertices in G{K). Let |$) = |^/''' (0)) be the same initial state as in Theorem 1. 
We choose the evolution time T so that for < t < T, the time-evolved state remains far from 
the vertices labeled {K, j) (for each j G {1, . . . , A^}), and thus far from the effect of truncating the 
paths. More precisely, we choose T = 0{L) so that, for times < t < T, the state \a^{t)) from 
Theorem 1 has no amplitude on vertices within a distance Nq = from the endpoints of the 
paths. For such times t we have 

(1 - P) H''\a^ (t)) = for all < r < Nq. 

This allows us to use the lemma with W = H = Hq\ |7(t)) = \a^{t)), and the bound 6 = 
0{L~^^^) from Theorem 1. The truncation lemma then says that, for < t < T, 



=o{L 



-1/4N 



(23) 



so, for < t < T, 



-iH 



(1) 



In other words, for small enough evolution times, the conclusion of Theorem 1 still holds if we 
replace the full Hamiltonian Hjj^ with the truncated Hamiltonian H^^^j^y 

We can also use the truncation lemma to extend Theorem 2 to the case where the infinite path 
has been truncated to a finite path. Let W = H = H^'^^ be the two-particle Hamiltonian (12) and 
define the truncated Hamiltonian H^'^'^ = PH^'^^P, where 



P 



\x,y){x,y\ 



x,y= 



with Ki,K2 = il{L). We take |$) = |^(0)) as in Theorem 2 and choose the evolution time 
T = 0{L) so that, for all times < t < T, the state \a{t)) has no amplitude on states where either 
particle is located within a distance Nq = n{L) from the endpoints of the truncated path. With 
these choices, and letting \'j{t)) = \a{t)) and 5 = 0{L^^/'^), we get that 



\m)-Ht)) =o{L 



for < t < T, which is the conclusion of Theorem 2 but now applied to the truncated Hamiltonian 
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E.2 Wave packet propagation on more complicated graphs 



The results of the previous section apply to wave packet scattering on finite graphs with long paths 
attached. The intuition behind these results is that, as long as the wave packets never get close to 
the ends of the paths, it does not matter whether we evolve using the Hamiltonian for the finite 
or the infinite graph. We now take this argument a step further. We expect that attaching another 
finite graph to the ends of the truncated paths will not substantially alter the time evolution as long 
as the wave packets never get close to the ends of the paths. 

We can use the truncation lemma a second time to make this intuition precise. For example, 
consider the case of a single-particle wave packet scattering on a graph G{K) of the form de- 
scribed in the previous section. Let |$) = \ijj^{0)) be the initial state of the particle as defined in 
Theorem 1, and choose the time T = 0{L) so that for all < t < T, \a^{t)) has no amplitude 
on vertices within a distance Nq = from the endpoints of the paths. We have already shown 
in the previous section that under these conditions the truncated Hamiltonian H^^jj^^ generates ap- 
proximately the same time evolution as the infinite Hamiltonian H^^^ (up to an error term that is 
0(L'^^^)). Let G* be a graph obtained from G{K) by attaching a finite graph to the endpoint 
vertices of G{K) (the vertices labeled (Z^, 1), . . . , {K, N)). 

Now apply the truncation lemma using W = H = H^^^j^^ and H = H^}, again letting P 

project onto G{K). With |7(t)) = \ci^{t)) and the bound 5 = 0{L^^/'^) from the previous section, 
the truncation lemma gives (for < t < T) 



G{K) 



l^^-(o)) 



H, 



(1) 



-1/4 



and hence 



O 



H. 



(1) 

G* 



^-1/4 



We see that the evolution of a wave packet for small times depends only on a portion of the graph 
(up to error terms bounded as above). This allows us to analyze a single-particle wave packet trans- 
mitting through a complicated graph consisting of subgraphs Gi{K) and G^iK) that overlap on 
long paths (e.g., the graph associated with the upper encoded qubit in Figure 13). For the purpose 
of analysis, we divide the total evolution time into intervals and evolve the state according to the 
Hamiltonian of a specified subgraph during each interval. Later we will see that this approach can 
even be used to approximate the time evolution of a system of more than one particle propagating 
through a graph. 

In Appendix E.3 we describe in full detail the finite graph implementing a single-qubit gate 
on one logical or mediator qubit. Then in Appendix E.4 we perform a similar calculation for the 
two-qubit C^^ gate. Finally, in Appendix E.5 we show how to combine these bounds to analyze a 
graph corresponding to an entire quantum circuit. 



E.3 Single-qubit gates 

Let us now apply the results of Appendix E.l to scattering on a graph implementing a single- 
qubit gate in our encoding. To approximate a single-qubit gate U on either a computational or the 
mediator qubit, we use a graph G{K) of the form shown in Figure 14, with 4 paths of length K 
extending outward from a subgraph G. We consider scattering at momentum k = — 7r/2 (in which 



33 



case U is either the identity or the Hadamard gate, up to a global phase) or A; = — 7r/4 (in which 
case U is either the phase gate, the identity gate, or the basis-changing gate, up to a global phase). 
We treat both choices for the momentum k in this section. By construction, the subgraph G has 
a 4 X 4 S-matrix of the form (5) at momentum k, with lower left submatrix equal to U. In this 
section we prove results about the single-particle evolution generated by the Hamiltonian H^^j^y 
We apply these results in Appendix E.5 where we analyze the full (n + l)-particle system. 

By choosing K to depend on the momentum k, we can compensate for differing particle speeds 
and have the scattering process occur over a fixed amount of time, ti = 3L/2. As discussed in 
Appendix D, we choose K(k) = 2M{k) + L where 



M(-7r/2) = L M(-7r/4) 



'3^2-2^ 
4 



L 



A logical input state a|0) + is encoded using single-particle states |Oin) and |lin) that only 
have support on the top left path and bottom left path, respectively. The encoded states differ 
depending on whether A; = — 7r/2or— 7r/4, and are defined as 

M{k)+L M(k)+L 
^ x=M{k)+l ^ x=M{k)+l 

Starting at t = from the superposition 

|V^(0)) = a|0in)+6|lin), 

the computation proceeds by evolving |^/'(0)) with the time-independent Hamiltonian H^^^^j^y cor- 
responding to a quantum walk on G{K). Here K = K{k) but we leave the argument implicit. 
After time ti, the state is 

\^K{ti)) = e-'''oU'\m)- 
The state is approximated by 

iV^out) = {Uma + Umh) |Oout) + {U^a + Uuh) |lout) (24) 

where |Oout) and |lout) are defined as 

M{k)+L M{k)+L 

lOout) = 6-2'*^^°^'=^ e''="|x,3) |lout) = e-^^*!™^'^^ ^ e'^'"|x,4). 

^ x=M(k)+l ^ x=M{k)+l 

Using the results of Appendix E.l, the error in approximating lipKih)) by |^out) goes to zero 
polynomially quickly as L grows: 

mKih))-\^,,,)\\ = OiL-'/'). (25) 

The effect of evolving the input state |^(0)) for time ti is depicted in Figure 12. 
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(K,\) (3,1) (2,1) (1,1)/^ (23) (3,3) (K,3) 

o • • • d o • • • o 
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(K,2) " (3,2) (2,2) (1,2)V JlA) (2,4) (3,4) " (KA) 



Figure 14: A graph G{K) used to perform a single-qubit gate on an encoded qubit. 
To prove equation (25), we write 

Wli^Kik)) - \^out)\\ < lll^i^(tl)) - \^o.{tl))\\ + lll^oo(tl)) - |V^out)|| . 

We then apply the bounds 

|||^oo(tl))-|^out)|| =0(^"'/') (26) 

and 

mKih))-\^^ih))\\ = OiL-'/'). ill) 

Equation (26) follows from Theorem 1. Applying the theorem with T = ti and M = M{k), we 
find 

a\a\ti)) - &|a'(ti))|| = C(L-^/^), (28) 

~[2tj smk\-M{k)-l 

J2 e*'^^(f/oo|x,3) + f/io|x,4)) 

x=-[2ti smfeJ-M(A:)-L 
M{k)+L 

J2 e^'" {Uoo\x,3) + t/io|x,4)) + 0{L-'/') 

x=M{k)+l 



M(k)+L 

e^'" (f/oi|x,3) + f/ii|a;,4)) + 0{L-^'^) 

x=M{k)+l 

where the 0{L^^/'^) error terms arise from approximating the upper and lower limits of the sum- 
mation. Comparing with equation (24), we see that (28) implies (26). 

The bound (27) follows from equation (23) (the truncation lemma with |$) = |V'(0)), H = 
Hff^^y H = H^^l^y W = H,6 = 0(L-i/4), iVo = M{k), and T = h). 

E.4 A two-qubit gate 

The C^^ gate is implemented using the graph shown in Figure 15. In this section we specify the 
logical input states, the logical output states, the distances X, Z, and W appearing in the figure, 
and the total evolution time. With these choices, we show that a CO gate is applied to the logical 
states at the end of the time evolution under the quantum walk Hamiltonian (up to error terms 



where 

|a°(ti)) 



and similarly 



ll^oo(tl)) 
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Figure 15: Graph G' used to implement the gate. The integers Z, X, and W are specified in 
equations (29), (30), and (31), respectively. 

that are (9(L"^/'')). The results of this section pertain to the two-particle Hamiltonian Hq) for the 
graph G' shown in Figure 15. 



The logical input states are 

1 

7z 



M{-l)+L 
x=M{-f)+l 



A/(-f)+L 



e-*4^|a;,2) 



for the computational qubit and 

A/{-f)+L 



x=M(-f)+l 



M(-f )+L 



- y 



2/=A/{-|)+l 
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j/=M{-f)+l 



for the mediator qubit. We define symmetrized (or antisymmetrized) logical input states for a, 6 e 

{0,1} as 

We choose the distances Z, X, and W from Figure 15 to be 

Z = AL (29) 
X = d2 + L + M (30) 

W = di + L + M (^-j^ (31) 
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where 



di 
d2 



4 



5L + 2di 5^ 
V2 2 



With these choices, a wave packet moving with speed a/2 travels a distance Z+2di + L = 5L + 2di 
in approximately the same time that a wave packet moving with speed 2 takes to travel a distance 

Z + 2d2 + L = 5L + 2d2, since 



t 



II 



5L + 2di 5L + 2d2 



= '-^ . 

V2 2 



We claim that the logical input states evolve into logical output states (defined below) with a 
phase of e*^ applied in the case where both particles are in the logical state 1. Specifically, 



e-^4''*n|00in) 



c,m 



e-^4''*n|01in)'='" 



e-^^G'*n|iOi„) 



c,m 



|llm)^'" 



|00o,,)^'" 
lOlout)^'" 
llOout)^'" 
lllout)^'" 



(32) 
(33) 
(34) 
(35) 



where, letting Qi = 2W + Z + 4 - M {-n/i) - L and Q2 = 2X + Z + 6 - M {-n/2) - L, 



lOout)" 



-ituV2 



x=Qi+l 



I lout) 

II \™ 
I J- out/ 



x,2) 



a-=Qi+l 



y=Q2+i 



y=Q2+i 



and lafeout)^'"' = Sym(|aout)16out)'"). 

Note that the input states are wave packets located a distance M(k) from the ends of the input 
paths on the left-hand side of the graph in Figure 15. Similarly, the output logical states are wave 
packets located a distance M{k) from the ends of the output paths on the right-hand side. 

The first three bounds (32), (33), and (34) are relatively easy to show, since in each case the 
two particles are supported on disconnected subgraphs and therefore do not interact. In each of 
these three cases we can simply analyze the propagation of the one-particle starting states through 
the graph. The symmetrized (or anti symmetrized) starting state then evolves into the symmetrized 
(or antisymmetrized) tensor product of the two output states. 

For example, with input state |00in)'^'™, the evolution of the particle with momentum — 7r/4 
occurs only on the top path and the evolution of the particle with momentum —n/2 occurs only 
on the bottom path. Starting from the initial state |Oin)'^ and evolving for time tu with the single- 
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particle Hamiltonian for the top path, we obtain the final state 

using the method of Appendix E.l. Similarly, starting from the initial state |Oin)™ and evolving for 
time til with the single-particle Hamiltonian for the bottom path of the graph we obtain the final 
state 

|0out)'" + O(L-i/^). 

Putting these bounds together we get the bound (32). 

In the case where the input state is |10in)'^'™ (or |01in)'^'™) the single-particle evolution for the 
particle with momentum — 7r/4 (or — 7r/2) is slightly more complicated, as in this case the particle 
moves through the momentum switches and the vertical path. The S-matrix of the momentum 
switch at the relevant momenta is given by equation (6). At momentum — 7r/4, the momentum 
switch has the same S-matrix as a path with 4 vertices (including the input and output vertices). 
At momentum — 7r/2, it has the same S-matrix as a path with 5 vertices (including input and 
output vertices). Note that our labeling of vertices on the output paths (in Figure 15) takes this into 
account. The first vertices on the output paths connected to the momentum switches are labeled 
(X + Z + 7, 3) and {W + Z + 2), respectively, reflecting the fact that a particle with momentum 
— 7r/4 has traveled W vertices on the input path, Z vertices through the middle segment, and has 
effectively traveled an additional 4 vertices inside the two switches. Similarly, a particle with 
momentum — 7r/2 effectively sees an additional 6 vertices from the two momentum switches. 

To get the bound (34) we have to analyze the single-particle evolution for the computational 
particle initialized in the state |lin)'^- We claim that, after time tn, the time-evolved state is 

|lout)^ + 0(L-l/^). 

It is easy to see why this should be the case in light of our discussion above: when scattering at 
momentum — 7r/4, the graph in Figure 15 is equivalent to one where each momentum switch is 
replaced by a path with 2 intemal vertices connecting the relevant input/output vertices. 

To make this precise, we use the method described in Appendix E.2 for analyzing scattering 
through sequences of overlapping graphs using the truncation lemma. Here we should choose 
subgraphs Gi and G2 of the graph G' in Figure 15 that overlap on the vertical path but where each 
subgraph contains only one of the momentum switches. A convenient choice is to take Gi to be the 
subgraph containing the top switch and the paths connected to it (the vertices (1, 2), . . . , {W^ 2), 
(1, 5), . . . , (Z, 5) and (X + Z + 7, 3), . . . , (2X + Z + 6, 3)). Similarly, choose G2 to be the bottom 
switch along with the three paths connected to it. The graphs Gi and G2 both contain the vertices 
(1, 5), . . . , (Z, 5) along the vertical path. Break up the total evolution time into two intervals [0, ta\ 
and [tajtii]. Choose t^, so that the wave packet, evolved for this time with H^^^, travels through 
the top switch and ends up a distance 6(L) from each switch, partway along the vertical path 
(up to terms bounded as 0{L^^/'^), as in Appendix E.l). With this choice, the single-particle 
evolution with the Hamiltonian for the full graph is approximated by the evolution with Hq^^ on 
this time interval (see Appendix E.2). At time t^, the particle is outgoing with respect to scattering 
from the graph Gi, but incoming with respect to 6*2- On the interval [ta, tn] the time evolution is 
approximated by evolving the state with Hq \ During this time interval the particle travels through 
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the bottom switch onto the final path, and at tu is a distance M(— 7r/4) from the endpoint of the 
output path. Both switches have the same S-matrix (at momentum — 7r/4) as a path of length 4, so 
this analysis gives the output state |10out)^'™ up to terms bounded as 0{L^^/'^), establishing (34). 
For the bound (33), we apply a similar analysis to the trajectory of the mediator particle. 

The case where the input state is |llin)'^'™ is more involved but proceeds similarly. In this case, 
to analyze the time evolution we divide the time interval [0, tu] into three segments [0, Ia], [tA, ^b], 
and [ts, tii\. For each of these three time intervals we choose a subgraph Ga, Gb, Gc of the graph 
G' in Figure 15 and we approximate the time evolution by evolving with the Hamiltonian on the 
associated subgraph. We then use the truncation lemma to show that, on each time interval, the 
evolution generated by the Hamiltonian for the appropriate subgraph approximates the evolution 
generated by the full Hamiltonian, with error 0{L^^^^). Up to these error terms, at times t = 0, 
t = tA,t = ts, and t = tu the time-evolved state 

• rr(2) , 

e-*^G'*|llin)"'°^ 

has both particles in square wave packet states, each with support only on L vertices of the graph, 
as depicted in Figure 16. 

We take Ga to be the subgraph obtained from G' by removing the vertices labeled ( [1.85L] , 5) 
,...,( [1.90L] , 5) in the vertical path. By removing this interval of consecutive vertices, we dis- 
connect the graph into two components where the initial state |llin)'^'™ has one particle in each 
component. This could be achieved by removing a single vertex, but instead we remove an interval 
of approximately 0.05L vertices to separate the components of Ga by more than the interaction 
range G (for sufficiently large L), simplifying our use of the truncation lemma. 

We choose = 3L/2. Consider the time evolution of the initial state |llin)^'™ with the two- 

(2) 

particle Hamiltonian Hqa ^^^^ "^^^ states |lin)'^ and |lin)™ are supported on disconnected 
components of the graph Ga, so we can analyze the time evolution of the state lllin)"^'™ under 

(2) 

Hq^ by analyzing two single-particle problems, using the results of Appendix E. 1 for each particle. 
During the interval [0, t^], each particle passes through one switch, ending up a distance Q{L) from 
the switch that it passed through and Q{L) from the vertices that have been removed, as shown in 
Figure 16(b) (with error at most 0{L^^^^)). Up to these error terms, the support of each particle 
remains at least Nq = Q{L) vertices from the endpoints of the graph, so we can apply the truncation 
lemma using H = H§), W = H = H^l, T = t^, and 6 = 0{L-^'^). Here P is the projector 

(2) (2) 

onto states where both particles are located at vertices of Ga- We have PHq, P = Hq^ since the 
number of vertices in the removed segment is greater than the interaction range G. Applying the 
truncation lemma gives 

= 0(L-^/^). 

We approximate the evolution on the interval [t^, t^] using the two-particle Hamiltonian H^j^, 
where is the vertical path (1, 5), . . . , (Z, 5). Using the result of Appendix E.l, weknow that (up 
to terms bounded as 0{L^^^'^)) the wave packets move with their respective speeds and acquire 
a phase of e*^ as they pass each other. We choose ts = 5L/2 so that during the evolution the 
wave packets have no support on vertices within a distance 6(L) from the endpoints of the vertical 
segment where the graph has been truncated (again up to terms bounded as C(L~^/*^)). Using H^^ 
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I 1 1 1 

di L M(-f) 

Figure 16: This picture illustrates the scattering process for two wave packets that are incident on 
the input paths as shown in figure (a) at time t = 0. Figure (b) shows the location of the two wave 
packets after a time = 3L/2 and figure (c) shows the wave packets after a time ts = tA + L. 
After the particles pass one another they acquire an overall phase of e*^. Figure (d) shows the final 
configuration of the wave packets after a total evolution time tu = {Z + 2di + L)/ \/2. 

(rather than Hq}) to evolve the state on this interval, we incur errors bounded as 0{L^^/'^) (using 
the truncation lemma with Nq = Q{L), W = H = H^l, H = H^^) , and 6 = 0{L~^/'^)). 

We choose Gc = Ga\ in the final interval [t^, ^ii] we evolve using the Hamiltonian again, 
and we use the truncation lemma as we did for the first interval. The initial state is approximated 
by two wave packets supported on disconnected sections of Ga and the evolution of this initial 
state reduces to two single-particle scattering problems. During the interval [t^, in], each particle 
passes through a second switch, and at time tn is a distance M{k) from the end of the appropriate 
output path. 

Our analysis shows that for the input state |llin)'^'™ the only effect of the interaction is to alter 
the global phase of the final state by a factor of e*^ relative to the case where no interaction is 
present, up to error terms bounded as This establishes equation (35). In Figure 16 we 

illustrate the movement of the two wave packets through the graph when the initial state is | lli„)'^'™. 

E.5 Block-by-block analysis of the full graph for a circuit 

In this section we discuss how blocks such as those shown in Figure 10 and Figure 11 act on 
encoded data. For each type of block, we first consider the time evolution generated by the Hamil- 
tonian i/biock for the block in isolation, i.e., with nothing connected on either side. The truncation 
lemma lets us use our results about i/biock to prove results about Hq^^\ the Hamiltonian for the 
full graph G where the input and output paths of the block are connected to other graphs. 
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E.5.1 Blocks applying single-qubit gates 



First consider a block of type I (as shown in Figure 10). The results of Appendix E.3 show that if 
the input and output paths of this block are not connected to anything, the block applies the correct 
logical single-qubit gates to each encoded qubit. Define single-particle logical input states l^^in)^ 
and output states |-2out)"' for each computational qubit (z E {0, 1} and j E {1, • • • ,n}) and the 
mediator qubit (j = m) as in Appendix E.3. 

Suppose that the state of the system at time t = is 

iV'(o)) = symi Y "d4)'---i4rkr'r 

\i'e{o,i}"+i 

Evolving the state |V^(0)) for time tj = 3L/2 using the Hamiltonian i/biock for the block we get 
li'ih)) = e~*^''"*^|V^(0)). Letting hi, /i2, . . . , /in+i be the single-particle Hamiltonians associated 
with the n -I- 1 encoded qubits (all supported on disconnected components of the graph), 

iV(ti)) = Sym I J2 age-''''''\ziy . . . e~^^"'^\zr'Xe-'^"+'''\zr'+^)'^ 
\i"e{o,i}"+i 

= Sym [ Y c^z{x'\Ui\z')\xlJ' . . . (a;"|?7jz")|<J"(a;"+^|Kn|z"+0|x^+i)- 
= \fi{U)) + 0{nL-"') 

where in the last line we have defined to be the logical output state with the appropriate 

unitaries applied. In the above we used the fact that 



-1/4N 



x=0 

as proved in Appendix E.3. Note that the error terms for each of the n + \ encoded qubits add 
linearly to give a total error 0{nL~^/^). 

Of course we are interested in the scenario where the input and output paths of this graph are 
connected to other graphs as specified in Appendix D. Let Hq'^^^ be the full Hamiltonian including 

these connections. We use the truncation lemma to show that e~*^G ^ *i \ip{0)) is approximated by 
|/i(ti)). We let P project onto the subspace of states where all n + 1 particles are located on vertices 
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of the graph within the block, H = Hq^^\ and W = H = i^biock-' For times < t < ti, we have 



\m) = Sym 



i*e{o,i}"+^ 



. . .e 



ihnt I n\n —ih„+it i n+l\ m 



W)) + W)) 



(36) 



where ||e(t))|| = 0{nL ^/^) and only has support on states where each particle is located 

at least a distance Nq = M(— 7r/4) from the endpoints of the input/output paths of the block (note 
thatM(-7r/4) < M(-7r/2)). In particular, 



(1 - P) (^H^G^^^^ ' |/i(t)) = for all < r < M (-7r/4) . 

Equation (36) again follows from the results of Appendix E.3 (using the results of that section to 
approximate e~*'*"'*|2;i„)'" for each w and t < ti). Applying the truncation lemma we then obtain 



o 



H, 



(n+l) 
G 



n 



^-1/4 



E.5.2 Blocks applying two-qubit gates 

Now consider a block of type II (as shown in Figure 11). Without loss of generality we assume 
that the computational qubit involved in the gate is the nth encoded qubit. We label the vertices 
on the input and output paths of the two qubits involved in the gate as in Figure 15. For the n — 1 
computational qubits not involved in the gate we label the states on each of the 2(n — 1) paths as 
\x,qy forj G {1, . . . ,n-l}, g e {0, l},andx G {1, . . . ,21V + Z + 4}. Herex = 1 is the leftmost 
vertex on each path. Define the input and output states for each of the — 1 computational qubits 
not involved in the gate as 



M(-|)+L M(-f)+L 
^^x=M(-f)+l y^:,=M(-l)+l 



J 



where Qi is defined in section Appendix E.4. Similarly, define the logical input and output states 
for the nth computational qubit and the mediator as in Appendix E.4. The state at time t = has 



'To ensure that PHq'^^^ P = i?biock we take L large enough that, for any two vertices i,j in the block that are 
associated with two different encoded qubits, there are no paths of length < C in the full graph G that connect i 
and j. This is true provided L > C (since any such path has to contain at least L vertices of another block that are 
part of a CO gate subgraph). Since C is a constant and L grows with the size of the computation, this requirement 
is automatically satisfied. Note that there are no paths of length < C that connect two vertices on different input (or 
output) paths for the same encoded qubit since the subgraphs we use to implement single-qubit gates have no paths of 
length < C between two input or two output vertices — see the beginning of Appendix D. 1. 
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the form 



1^(0)) = sym Yl ■ ■ A^-Y-'KrKn'T 



, ^e{o,l}"^ 



\ze{o,i}"+i / 



V2 

where in the last line we have written the state in terms of the symmetrized (or anti symmetrized) 
logical states |00in)"'™, |01in)"''^, |10in)"''", 111^)"'"" as defined in Appendix E.4. The prefactor of 
1 / -\/2 arises from the fact that these logical states for the two encoded qubits n and m are already 
symmetrized (or antisymmetrized). 

Evolving this state for time tu using the Hamiltonian of the block (without anything connected 
on either side) gives |^/'(tn)) = e~*^'''°=''*"|?/'(0)). Let hi, /i2, • • • , hn-i be the single-particle Hamil- 
tonians for the components of the graph corresponding to the n — 1 encoded qubits that are not 
involved in the gate. Let /i„ m be the two-particle Hamiltonian for the two encoded qubits on which 
the gate acts. Then 



mtii)) 



^ \2e{o,i}"+i 



age 



-ihitji I ^1 \ 1 
^in/ 



. . . e 



n n+l\7 
^in / 



Using the results of Appendix E.4 this is 



\ze{o,i}"+^ 



= |/€(t„)) + 0(nL-i/^) 

where is the encoded logical state with the unitary applied. 

Now connecting this block to other blocks on either side, and letting Hq~^^^ be the full Hamil- 
tonian, we can apply the truncation lemma just as we did in the previous section for blocks of type 
I. This gives 



e-^r^^*n|^(0))-|«:(tn)) 



E.5.3 Piecing blocks together 

We are now ready to show how to apply two unitaries in series, by concatenating blocks of type 
I or II as shown in Figure 10 and Figure 11, respectively. We concatenate two blocks B and B' 
by removing some of the vertices on the input paths of B' and then connecting the output paths 
of B to the resulting graph (so that the two blocks overlap on the removed vertices). We remove 
the leftmost 2M{k) + L vertices from each input path of B' (this is a different number of vertices 
depending on whether the paths are associated with a mediator qubit or a computational qubit) and 
then connect every output path of B to the corresponding shortened input path in B' . Note that 
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if B' is of type I, we remove all vertices from each input path, attaching the output path of B to 
inputs of the individual gate subgraphs of B' . 

Assuming that our wave packet initially starts as the correct input for a single block (i.e., a 
symmetrized logical state with one particle on each pair of qubit paths), we want to show the wave 
packet first moves through B, undergoing the appropriate unitary, and then moves through B' , 
undergoing the second unitary, with small total error. Define ti and t2 to be tj or tn depending 
on the block type of B and B' , respectively, and let the initial state of the system be IV'(O)), a 
symmetrized logical state as described previously. From the results of the previous two sections, 
the state of the system after time ti is |outi) + ||r2L^^/^). Here |outi) is a symmetrized 

logical state in which we have applied the unitaries of B to the initial logical state. Furthermore, 
each individual wave packet is located a distance M{k) from the end of B, and has length L. Since 
B and B' overlap on 2M{k) + L vertices, we see that each wave packet is located a distance M{k) 
inside of B' , so the output state for B corresponds with an input state for B' up to an irrelevant 
global phase. After time ti + t2, the wave function is then jouts) + UnL-^/^), where 

|out2) is a logical state in which the unitaries of B' and B have been applied in sequence. Now 
generalizing to the case of g blocks in series, we see that the total error in the final output state is 
0{g\\HQ'^^^\\nL~^/^), as claimed in equation (19). 

F Single-particle wave packet scattering on infinite graphs 

In this section we prove Theorem 1 . The proof is based on (and follows closely) the calculation 
from the appendix of reference ((5). 

Recall from (4) that the scattering eigenstates of Hq'' have the form 



for each k G (— tt, 0). 

Before delving into the proof, we first establish that the state \ci^{t)) is approximately normal- 
ized. This state is not normalized at all times t. However, {a^ {t)\a^ [t]) = 1 + 0{L^^), as we now 
show: 



(x, g|sc,(fc)) = e-''="5g,- + e*'="5,,(A:) 



1 



oo 



{a={t)\a^{t)) 



L 



\e-"'''R{x - [2tsinA;J) + Sjj{k)e"''' R{-x - [2tsinA;J) 



x=l 




q^j x=l 



1 



oo 



L 



^ [R{x - [2t sin k\ ) + R{-x - [2t sin k\ )] 



x=l 



+ 




x=l 




x=l 
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where we have used unitarity of S in the second step. When it is nonzero, the second term can be 
written as 



1 J2 (e-2^*=-5^(fc) + e^'''^S,,{k)) 



x=l 



where b is the maximum positive integer such that {—b, b} C {M + 1 + [2t sin k\ , . . . , M + L + 
[2t sin k\ }. Performing the sums, we get 



x=l 



1 

L 



^-2ikb _ 



< 



LI sin A; I 



Thus we have {a^ {t)\a^ {t)) = 1 + 0{L-^) 
Proof of Theorem 1 . Define 



and 



N 



q=l 



where we take e = Observe that is a projection (i.e., = IIJ, as can be shown using 

the delta- function normalization of the scattering states. Thus we can write 



iV^^'(t)) = \w^it)) + \v\t)) 



where 



\w'it)) = n,|^^'(t)) 



AT 



^g-2.*cos(fe+0) J2 \sc,{k + 0))(sC,(fc + <P)\ij\0)) 
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q=l 



and {w^ {t)\v^ {t)) = 0. Now 

(sc,(A: + 0)|^^-(O)) = 4r E (^^"^'^^^^ + e-'^"'^'^^^S*Jk + 

^ ^ x=AI+l 



M+L 



SO 



N 



\wKt)) = WA{t)) + 

5=1 
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where 
with 



f(^^ - -If E ^^'^ 

x=M+l 

M+L 

We will see that iV'^Xi)) ~ k^W) ~ 

We have 

K(t)K(t)) = y_^- 1/(0)1 =zy_,^^Ii?(i^' 



but 
and 



1 r rf0sin2(iL0) _ ^ 



1 
L 



/ r ^siii^(|L0) _ 2 r d0 sin^(|L0) 

;U ^J-j2^ sin2(l0) 2^ sin2(l0) 

2 rrf07r2 



2 r d(f)7r^ 



Therefore 
Similarly, 



TT 
< 

- Le 



i>K(t)K(t))>i--^. 
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and, using the unitarity of S, 



N 



q=l 



d(P sm\^L{2k + (P)) 
L 2^ sin2(i(2A; + 0)) 



< 



Now |sm(A; + 0/2) — sin A;| < |0|/2 (by the mean value theorem), so 



sin^ (^"^2) ~ (^N^^^l ~ 



Since e = < \sink\ we then have 



Hence 



^ 1 



27r sin^ k 



4e 



TrLsin /c 



K(t)|^.^(t))>K(t)K(t))-2 



N 



TT 



9=1 



9=1 



q=l 



> 1 4'i/ 5—, 

-Le V TrLsin'^A; 



so 



Le V TrLsm A; 



since {v^ {t)\v^ (t)) + {w^ {t)\w^ {t)) = 1. Thus 



||IV'^'W)-KW)lh 



N 



9=1 



< I ^ + 4 / 
~ Le V ttL sin^ /c 



eiV 



TrLsin /c 
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With our choice e = -1^^, we have — = ^(-^ ^^^)- We now show that 



(38) 



Letting 

N oo 

P = ^^\x,q){x,q\ 

q=l x=l 

be the projector onto the semi-infinite paths, to show equation (38) we use the bounds 



[i~p)\w'At))\\ = o 



\ogL\ 



and 



\P\w\it))-\a^it))\\ = 0{L-'/'). 



(39) 



(40) 



Equation (38) follows from these bounds since (1 — P)\a^{t)) = 0. 
To get equation (39), write 



{wm^-p\<it)) 



ID, 
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Using Lemma 1, there is a constant A such that \{scj{k + 0)|1 — P\scj{k + 0))| < A^m, so 



' — e J —e 



^1/(0)1) . 



(41) 



Now 



sm 



(f>L 



smf 



< 



< 



Jo 27rVL 



sm 



<j,L 



sm| 



for any d E (0, e] 



+ 2 / r- 

TT 7d 27rVL' 



TT 



dVL , logj 



+ 



(42) 



where in the last line we used the fact that e < 1. Setting d = 6(1/L), using equation (41), and 
taking the square root of both sides, we get equation (39). 
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We now prove the bound (40). Noting that 

1 



we write 



{x,q\a^{t)) = e 



- 2it cos k 



ikx 



-i4>{x— [2t sin 



27r 



'=J)/(0) 



ikx 



^g-^<^i(-x-L2^sinfcJ)J(^0^j j _ 



On the other hand, 

(x,g|w^(t)) 



^^-2Ucos(k+cj,)f^^^ ^^^^ik+<p)x^^^ ^ e*(^'+*)-5,,(A; + 0)) 



Using equations (43) and (44) we can write 

PK(t)) = |a^-(t))+X:icl'(t)) 

where P|c^(t)) = |4(t)) and 



i=l 



{x,q 
{x,q 
{x,q 
{x,q 
{x,q 
{x,q 
{x,q 



4{t)) = 6, 



-2it cos k ^— ikx 



4W) = s,,{k) 



-2it cos k ikx 



-2it cos k^—ikx 



^gj0a'y(^0) ^g2ii0sinfc _ gj0[2tsinfcj^ 

+ 



27r 



g-j0zj^0^g2jt(^ sinfc 



= -5,,(fc)e 



-2it cos k Akx 



7T P — ^ 
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g^i<l>^ y(^(^^g2it0 sinfc 



4{t)) = s. 



-ikx 



27r 



-i(/>x'y^0^ 2jt cos(fe+</') g— 2it cos fc+2it()!> sin fc^ 



4W) = 



4W) 



,jfcz 



27r 



27r 



-2jtcos(fc+0) — 2ii cos fc+2it(/) sin fc'^ 
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We now bound the norm of each of these states: 

N oo 



q=l x=l 

N oo 



6gje 



-2it cos k^—ikx 



^2itif>s'mk i(/> [2t sin fc J ^ 



2n 



q=l x= — oo 



6gje- 



-2it cos k ^—ikx 



^^g-i(/>xy^0^ ^g2it0sinfc _ gi(^[2t sin fcj j 



2 I 2ii0sinfc _ i^[2tsinki 1 2 



— e 



</ — |/(0)r(2t0sinA;- [2tsinA;J 



2 i,2 



where we have used the facts that le*"* — 1 p < for s G M and \2t sin k — [2t sin A;J | < 1. In the 
above we made the following replacement under the integral: 

oo 

^i{4>-4')x ^ 27f 6{(f) - 0) for 0, G {-71, tt). 

x=— oo 

We use this repeatedly in the following calculations. Continuing, we get 



(ci(t)ici(t)) < 



1 r d(t)sm.WL<i)) 



L ./_^ 2-n sin^(|( 



1 

< - 



7r2 
< — 
- L 



using the fact that sin^(0/2) > ^^/tt^ for G [-tt, tt]. Similarly we bound < vr^/L. 

Using equation (37) we get 



+ 



Zti 



TT 

< 

- Le 
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and similarly for {c'^{t)\c^^{t)) . Next, we have 

{4{t)\4{t)) < y ^ 1/(0)1' |e-2^*^°^('=+*) - e-^it CO. k+2it<t>.mk^'' 

< f — \f {(I)) f {2t COS {k + (j)) -2tcosk + 2t(j) sink f 
J-e 27r 

= / — 1/(0) I' (2t COS A; (cos 0-1) + 2t sin A; (0- sin 0))^ 

< £^1/(0)^4^^ 



4^2 



L J _^ 271 sin^d^ 



and we have the same bound for |(4(t)). Finally, 

N 



(4(t)i4(t)) < r ^ 1/(0)1' |5,,(^ + 0) - s,,ik)t 

^ q=l 

Now, for each g e {1, . . . , A?^}, 

|5,,(fc + 0)-5,,(/c)|<r|0| 



where the Lipschitz constant 

d 



r = max max 

q,je{l,...,N} fc'e[-7r,7r] 



is well defined since each matrix element Sqj{k') is a bounded rational function of e*'^' (see Section 
3 of (54)). Hence 



{4(t)m)< £g 1/(0)1' TvrV 



NT^ f rf0sin2(lL0) 
_,2^^h?(|0)" 
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Now using the bounds on the norms of each of these states we get 



using the choice e = and the fact that t = 0{L). □ 



G Two-particle wave packet scattering on an infinite path 

In this section we prove Theorem 2. The main proof appears in Appendix G.l, relying on several 
technical lemmas proved in Appendix G.2. The proof follows the method used in the single- 
particle case, which is based on the calculation from the appendix of reference (6). 

Recall from (16) that for each pi E (— tt, tt) and p2 G (0, tt) there is an eigenstate |sc(pi; P2))± 
of H^^^ of the form 

_ipi(£±2) f e-*P2(^-?/) ± e»''±(pi'P2)e^P2(^-y) [f ^ - y < -C 
{x,y\sc{pi;p2))± = ^ ^ < e-*P2(^-?/)e^^±(Pi'P2) ± ev^^^-v) lix-y>C (45) 

[/(Pi,P2,a; - y) ± f{pi,P2,y- x) if \x-y\<C 

where 

e''^^P''P''> =T{p,,p2)±Ripi,p2), 

C is the range of the interaction, T and R are the transmission and reflection coefficients of the 
interaction at the chosen momentum, / describes the amplitudes of the scattering state within the 
interaction range, and the ± depends on the type of particle (+ for bosons, — for fermions). The 
state |sc(pi;p2))± satisfies 

i?(^)|sc(pi;p2))± = 4cosy cosp2|sc(pi;p2))± 
and is delta-function normalized as 

±{sc{p[;p'^)\sc{pi;p2))± = ^7r^S{pi - p[)S{p2 - P2). (46) 



G.l Proof of Theorem 2 

Let 

in 



= / / A 2 ^ |sc(pi + 0i;P2 + 02))±±(sc(pi+0i;p2 + 02)| 



with D, = [-e,e] x [-e,e], pi = n /2 - n /A = n/A, and ^2 = (7r/2 + 7r/4)/2 = Svr/S. By the 
delta-function normalization of the scattering states (equation (46)), 11^ is a projection. Thus we 
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can wnte 

\^it)) = e-^''''''\m) = mt)) + mt)) 

where 

|^i(t)) = n,|^(t)) 

±(sc(pi + 0i;p2 + 02)1^(0)) 

and |V'2(^)) is orthogonal to We take e = a/ for some constant a. Using equation (45) 

we get 

\Mt)) = \^A{t))±\Mt)) 

where 



5(01, 02, f ) |sc(f + 0i; f + 02))± (48) 



with 

-(M+l) Af+L 

A(0i,02) = ^ e*0i^e^<^^(^-^) (49) 

a:=-(Af+L) y=M+l 



x=-{M+L) y=M+l 

Again using the delta-function normalization of the scattering states, we get 

{Mt)\Mt)) = ll/-^\B{<p., 

167r2 
< 

by Lemma 5 (provided e < Svr/S, which holds for L sufficiently large). Similarly, 

i>{Mt)\Mt)) 

. 2 1^(01' 02)1 
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(from the first two facts in Lemma 5) and therefore 

(MtMiit)) = {Mt)\Mt)) + {Mt)\Mt)) + {Mt)\Mt)) + {Mt)\Mt)) 
>i-^-2\{MtMBm 

Le 

> 1 _ g _ 2 |(^ij(t)|^B(t))|^ 

127r 

> 1 



Hence 



Le 

12tt 



smce 



Thus 



Le 



II IV^W)-I^^W)ll = ll l^i.W) + l^2(t))|| 

<ll I^B(t))|| + || |V'2(t))|| 



Now 



i*(«)> - ia(())ii < II \m) - \Mtm + II \Mt)} - Hm 



An 1271 ,, , , , , , , , , , ,, 
-L^ + VTr + " I^AW)-|a(t))|| 
= 0(L-V4) + || |^^(^))_|^(^))|| 

using our choice e = a/\/L. To complete the proof, we now show that the second term in this 
expression is bounded by 0{L^^/'^). 

Lemma 3. With and |a(t)) defined through equations (47) and (20), respectively, and with 

t < cqL (for some constant Cq), \\\'ipA{'t)) — |"(^))|| = 0{L~^/'^). 

Proof. To simplify matters, note that for x ^ {x,y\a{t)) = ±{y, x\a(^t)) and {x,y\ipAit)) = 
±{y, xl^Ait))- Taking C to be the maximum range of the interaction in our Hamiltonian, we have 

II - Hm < 2 \\Pi\Mt)) - PiHm + ii^2i^A(t))ii + iiP2i«(t))ii , 

where 

^1 = k,y)(a;,l/| P2= l^'^)(^'^l- (50) 

y-^>c \x-y\<C 
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Now, fory — x>C, 



{x,y\iJA{t)) 



2 J 



d(f)id(f)2 
47r2 



V2 



g-i7rx/2gi7r?//4g-i</>i(^)gi</>2(y-x) 
^^inx/4^-i7ry/2^-i<t>i{^)^-i<jy2{y-x)^ie±{^+<t>u^+^2) 



From Lemma 6, for x < y, the state takes the form 



g-i7rx/2gi7ry/4 



where Z^^r = [— vr, vr] x [— tt, tt]. Using these expressions for \tpAit)) and we now write 

PilV^^(t)) - Pi|a(t)) = ±|ei(t)) + |e2(t)) ± |/i(t)) + l/sW) ± \giit)) + |^?2(t)) ± m)) 
where each term in the above equation is supported only on states \x, y) such that y — x > C, and 
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(for y-x>C) 
{x,y\ei{t)) 



{x,y\e2{t)) 



{x,ym)) 



{x,y\f2{t)) 



{x,y\g2{t)) 



{x,y\h{t)) 



-t<pi 



_}_^-itV2^-iTrx/2^my/4: 

V2 



1^ 



^(01,02) 



/■ /■ d(j)id(j)2 



_^_^-itV2^iTrx/4:^-iTTy/2 
V2 



47r2 



-^(01,02 



g it\/2^-iTTx/2^iny/4 



d(f)id(f)2 



47r2 



-^(01,02) 



(x,t/|^i(i)) = -=e--/^e 



inx/A —iiTy/2 



V2 



,-it4cos(f + ^)cos(^+02) _ -it(V2+V2(^-<P2)-2(^+4>2)) 



_}_^-iTrx/2^iTTy/4 



V2 



^ g j7ra;/4g-i7rj//2 



^A{<p„<P2)e-M'^)e-'M^: 



V2 



g-it4cos(f + ±L)cos(^+</,2) _ ^-it(V2+V2{!!^~4>2)-2{^+4>2)) 



We now proceed to bound the norm of each of these states. We repeatedly make the following 
replacement inside the integrals (here 0i, 02, 0i, 02 e (— tt, tt)): 

oo 

^ gix(i(<^l-^l)-(<^2-^2))gi2/(i(<^l-^l) + (<^2-02)) ^ 47^2^(0, _ 0i)5(02 - 02). 
a;,2/=— oo 
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Using this formula we get 



(ei(t)|ei(t)) = {eiit)\x,y){x,y\eiit)) 

y-x>C 
oo oo 

a;=— oo y=—oo 



1 ff d(t)id(t)2 , . 



72 ' 2 



d(t)id(t)2 . , .,2 

1^(01,02)1 



47r2 



-<^i(-W + [^J)g-2i02(-LtJ-[ 



%/2 



Now use the fact that |e — 1 1 < for c G M to get 



(ei(t)lei(t)) < \ [^^^) 1^(01,02)1^ -01 (^-^ + 71 + W 



.2 1 -t - ^ + L^J + 



< 4 



iA(0i,02)r(0?+40^) 



using the Cauchy-Schwarz inequality and the fact that |t — t/y/2 — [t\ — \t/\/2\ | < 2. So 



(ei(t)|ei(t))<4 



d(f)id(f)2 

47]-2 



+ 



yi»y2 
47r2 



|A(0i,02)r(0? + 40^) 



< 4 ■ Svr^— + 20e2 
Le 

where we have used Lemma 5 and the fact that 0^ + 40| < 5e^ on D^. Similarly, 



Now 



(e2(t)|e2(t)) < 



1 



SOtt^ 



+ 20e' 



{m\m)<- 11 

< 

- Le 



^"'^^ 1^(01,02)^ 
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by Lemma 5, and similarly 



Moving on to the next term, 



(/2WI/2(t)) < 



27r 



{9i{t)\giit)) < - 



1 

< - 
- 2 



47r2 



i^(0i,02)r 



- jt4 cos( f + ^ ) cos( +02 ) 



-it(V2+V2(^-<^2)-2(^+<^2)) 



47r2 



|A(</)i,</)2)|'t2(4cos(- + :|)cos 



using |e — 1| < for c G M. Now 



fir 01 \ /Stt 

f -K (hi \ f -K (hi 

= 2 cos f - + ^ + 02 j + 2 cos f -- + ^ - 02 



-2 sin (^^ + 02^ + V2cos(^y 



+ V2sin ( ^ 



so 



7r 0i\ /37r 
4 cos I — + — I cos 



< 



%/2 



cos 



- 1 



)-v^-v/2(^-02) 
\/2 ( sin 



+ 



2 I sin I Y + </;2 



+ 



+ 



^2 I + 2 ( y + </^2 
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using |cosx — 1| < and |sinx — x| < for x G M. Plugging this into equation (51), we get 



{gi{t)Mt))<- 



2 



16t2 rr rf0^(i02sin'(f[^ + 02])sin2(f[-^ + 02]) 



JJd^ 4vr2 sin2(l[f + 02])sin2(i[-f + 02]) 



1 , , \ , / 01 



< 



y + 02j 

where we used the Cauchy-Schwarz inequality in the second line and equation (55) in the last line. 
Changing coordinates to 

ai = 01 + Y ^2 = Y - 02 

and realizing that |q;i|, |a2| < 3e/2 for (0i, 02) G D^, we see that 

{gi[t)\gi(t)) < -— — . 271 X — r^T^^ "i + "2) 

^ i-3e/2 27r y_3,/2 SvT SIH (ifti) SIR (^^2) 

32t2 f^'/^ dai f^'/'^ da2sm'^{Uai)sm^ULa2) 4 



-^^ 7-3£/2 27r 7_3,/2 27r sin (2ai) sin {^02 
^ 32^2 /'3^/2 dai sin^dLai) ^ 
^ 7-3C/2 27r sin^(iai) ^ 
32^2 y3^/2 da, 7r2 

367rt2e3 



L y_3e/2 27r a2"i 



with the same bound on {g2{t)\g2{t)) . 
Finally, 



1 



Recall that e^^±(Pi'P2) = T{pi,p2)±R{pi , P2) is obtained by solving for the effective single-particle 
S-matrix for the Hamiltonian (14). For pi near tt/A we divide this Hamiltonian by 2 cos(pi/2) to 
put it in the form considered in (34), where the potential term is now V(|r|)/(2 cos(pi/2)). The 
entries T(pi,p2) and R{pi,p2) of this S-matrix are bounded rational functions of 2; = e*^^ and 
(2 cos(pi/2))^^ (54), so they are differentiable as a function of pi and p2 on some neighborhood 
U of (7r/4, 37r/8) (and have bounded partial derivatives on this neighborhood). 

For e small enough that Z^^ = (7i"/4, 37r/8) C U we get, using the mean value theorem and the 
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fact that = 0±(7r/4,37r/8), 



gi6»±(|+01,^+02) _ gj' 



< \J(tA + 0i max I Ve^'^i | for (/.g) G 

< er 



for some constant F (independent of L). Therefore 

< -rV. 

- 2 

Putting these bounds together, we get 

\\PMA{t)) - PMm < ll|ei(t))|| + |||e2(t))|| + + |||/2(t))|| 

+ \\\9im + \\\92m\ + mm 

Letting e = a/ vT and t < CqL we get 

II PiI^aW) - Pi|«(t))|| = OiL-'/^). (52) 

Since P2|ct(t)) has support on at most ACL basis states \x,y), and since \ {x,y\P2\a{t))\'^ = 
we get 

||P2|«(t))|| = 0(i^-'/'). (53) 



Furthermore, Lemma 7 says that 

/ 

II pj,i,Jf\\\\ ^ n I : 



P2|^A(t))||=0(^). (54) 



Using equations (52), (53), and (54), we find 

IIIV^A(t)) - |«(t))|| < 2 IIPiIV^aW) -Pl|«(t))|| + ||P2|«(t))|| + ||P2|^A(t))|| 

= 0(L-i/"), 

which completes the proof. □ 
G.2 Technical lemmas 

In this section we prove three lemmas that are used in the proof of Theorem 2. 
Lemma 4. Let \<y{t)) be defined as in Theorem 2. Then 

{a{t)\a{t)) = 1 + 0{L-^). 
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Proof. Define 

x<y 

Note that, since {x, y\oi{t)) — ±(y, x\a{t)) for x ^ y. 



{a{t)\a{t)) ^ 2{a{t)\Ii\a{t)) - {a{t)\x,x){x,x\a{t)) 

x=—oo 

= 2{a{t)\U\a{t)) + 0{L-^) 

where the last line follows since | (x , x | a (i) ) p is nonzero for at most L values of x and | (x , x | a (i) ) | ' 
C(L-2). We now show that 



{a{t)\U\a{t))^- + 0{L-'). 



Note that 



{a(t)\U\a(t)) = ^Y.i ^(^' ^' ^) + ^(^' ^' ^) 



3^7r 3i7r 



± e'^e— ^e"— y, t)F{y, x, t) 



Now 



and 



J2F{y,x,t) = Y,F{x,y,t) 

x<y y<x 



^ ^ / 00 00 00 

a;<2/ \x=— 00 y=—oo x=—oo / 



We now establish the bound 



_1_ J2 e'^^e-''^yF{x, y, t)F{y, x, t) 



x<y 



to complete the proof. To get this bound, note that both F{x, y,t) = 1 and F{y, x,t) = 1 if and 
only ifx,y G B where 



B = {-M-L + 2[tJ,...,-M-l + 2[tJ}n|M + l- 2 



V2J' 



,M + L-2 



LV2. 
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Observe that 



for some j,l G Z with I < L. So 



B = {j,j + l,...,j + l} 



1 
2Z2 



J2e'^^e-'^yF{x,y,t)F{y,x,t) 

x<y 



e 4 4 

x,y&B, x<y 

j+i y 



E 



< 



2L2 
1 

2Z2 
1 

y=3 
/ + 1 2 

2L2 |e3^f - l| 

= Oil-') 

since I < L. 

Lemma 5. Let k G (— tt, 0) U (0, tt) anJ < e < min {vr — \k\, \k\}. Let 

De = [-e, e] X [-e, e] 
= [-n,7r] X [-7r,7r] . 



e3^f _ 1 



Then 



47]-2 



1^(01, 02)P 



4772 



■|i?(0i,02,A;)|^< 



where ^2) ^inJ -B(0i, 02? ^) ^ire ^/ven by equation (49). 
Proof. Using equation (49) we get 

-(M+l) M+L 



^i^{x+y-{x+y)) ^i<j>2{x-y-{x-y)) 



Now 



a;,i=-(Af+L) y,y=M+l 



<MH Mx-y~x+y) _ c 
2^ C — Ux-y,x-y, 
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so (suppressing the limits of summation for readability) 



x,x y,y 



x,x y.y 



which proves the first part. 

By performing the sums in equation (49) we get 

Letting ai = </)i/2 + (f)2 and 0^2 = (^i/2 — 02, we see that |q;i| < 37r/2, |q!2| < 37r/2, and 
a? + ai > 5eV2 whenever (0i,02) e \ A- Defining D^^i2 = [-37r/2, 37r/2]2 we get 
(ai, 02) e -D3,r/2 \ A whenever (0i, ^2) ^ D^. Hence 



2^ 1 [[ ciaida2sin2(|LQ;i)sin2(lLa2) 
- L 271 sin^dai) J ^ 27r sin2(la2) 



<-( r 1 sin^(i-^Qi) \ / /"^ c^a2 1 
- L "271 sin^dai) ) ^^h?(I^ 

_8 f r da2 1 /"^ rfa2 1 

L yi, 27r sin2(ia2) A Svr sin2(ia2) 



2ti al ' " 27r 



47r 



which proves the second inequality (in the next-to-last line we have used the fact that sin(x/2) > 
x/7rforx e (0,7r) andsin^(x/2) > 1/2 for x G (7r,37r/2)). 
Now 

|S(01,02,fc)P = |A(0i,-02-2fc)|2 

1 1 1 

- Z^sin2(i[f + 02 + 2k]) sin2(l[-f + 02 + 2A;]) ' 
If (01, 02) G -D, then -3e/4 < |±0i/4 + 02/2 + A;| < |A;| +3e/4. Noting that e is chosen such 



63 



that < e < min{7r — \k\, \k\}, we get 



< 



4 2 



< 77 - 



so 



15(01, 02, A:) P< 



< 



1 1 



and 



|5(0i,02,A;)| < — (2e) 



47r2 



47r2 ' ' V^^eV 



L2e2 



Lemma 6. Le? a^yit) be as in Theorem 2. For x < y, 



^(01,02) 



47r2 



47r2 



^(01,02) 



□ 



Proof. The lemma follows from (21) and the fact that, for any two numbers 71, 72 such that 71 + 

72,71 - 72 e z, 



47r2 



A(0i,02)e 



j7i</'l+2i72 02 — ) L 



J if (-71 -72,-71 + 72) e ^5 



otherwise 



where S = {-M - L, . . . , -M - 1} x {M + 1, . . . , M + L}. To establish this formula, observe 
that 



-M-l M+L 



^ ^ U(0i,02)e^^^<^^+2n202 



L ^ ^ J ^ 2ti ^■ 



'X— 272 • 
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Here we have performed the integral over 02 using the fact that 272 is an integer. We then have 

-M-l M+L 



L ^ 2^ y 27r ^y,— 271 



x=-M-L y=M+l 
'M~l M+L 



x=-M-L y=M+l 



as claimed. 



□ 



Lemma 7. Lef Z^e a* m equation (47) with e = (for some constant a), and let P2 be as 

in equation (50). Then 



Proof. First note that 

KMtmiMm < 

Now 



d(f)id(f)i f d(f)2d(f)2 
47r2 



47^2 



A(0i,02)M(0i,02) 



±(sc(f + 0i; f + 02)|P2|sc(f + 0i; f + 02))± 



±(sc(f + 0i; f + 02)|P2|sc(f + 0i; f + 02))± = ±(^A(f + 0i; f + 02)| J|^(f + 0i; f + 02))± 
where 

m-, + 0i; f + 02))± = ^ (|^(f + 0i; f + 02)) ± |^(f + 0i; -f - 02))) , 

P2)) is defined in equation (15), and 

J=^e<^-^Of J2 |r)(r| + 5^e<*-^0f ^ |r)(r| 

seven reven, |r|<C s odd r odd, \r\<C 

= 27r5(0i-0i) ^ |r)(r|. 

\r\<C 

Define 

^(01, 02, 02) = Yl + f + 02)|r)(r|V^(f + 0i; f + 02))±. 

|r|<C 

The function g satisfies |5'(0i, 02, 02)1 < 4CA(0i)^, where for each 0i, A(0i) is the constant 
obtained by applying Lemma 1 to the single-particle Hamiltonian 

-1 



i/«+(2cos(f + f))-^5^V(|r|)|r)(r|. 
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Let p = AC max[_f_e] A(0i)^. Then 

\{i^A{t)\P2\^A{t))\< 



d(j)id(f)2 p 



271 



yl((/)l,02)*^(01,02 



47r2 L 



sm(|(i0i + 02)) sm(f - 02)) 



sin(i(i0i + 02)) sin( 



01 - 02)) 



2711 



sm(f (i0i + 02)) sm(f (i0i - 02)) 



sin(2(20i + 02)) sin(2(20i - 02)) 
d(j)id(j)2 



A7T^L 



sin(f (i0i + 02)) sin(f (i0i - 02)) 



sm 



1 f 1 



,2V2^1 +02)) Sin(2(201 - 02)) 



AtiL 



sm 



-(- 

. 2 V2 



01+02)) 



sm 



-(- 

2V2 



01 + 02)) 



+ 



sin(f(|0i -02)) 
sin(i(i0i - 02)) 



Now for each 0i G [— e, e], 



sin(f(i0i±02)) 



sm 



2(2 



01 ±02)) 



< 



< 



I' dd 



27iL 



sin(f02) 



sm 



27iL 



sin (1^02 



i02) 

2 



27tL 



sm 



sin(^02) 



sm 



so 



A7r^L 



sm(f (|0i + 02))sin(f (i 



01 - 02)) 



sin(i('i0i + 02)) sin( 



>2V2 



|(|V1 



01 - 02)) 



(56) 



Letting ai 



>! + 02 and 02 



)i - 02, we get 



1 
L 



47r2 



sin(f (i0i + 02)) sin(f (|0i - 02)) 



sin(2(20i + 02)) sin(2(20i - 02)) 



< 



< 



dai 
~2^ 



sin (ftti) 



sin 



dao 



2ti 

2*= 



sin (f 0:2) 



sin (|a2) 



' d^ , logi^ 



+ 



TT 



for G (0 



-1 

' 2 J 



where in the last line we used equation (42). Setting d = 0(1/ L) and using this bound in (56) 
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gives the desired result. 



□ 



H Truncation lemma 

To prove Lemma 2 we use the following two propositions: 

Proposition 1. Let H be a Hamiltonian acting on a Hilbert space Ti, and let 1$) G Ti be a 
normalized state. Let K be a subspace of Ti such that there exists an Nq E N so that for all 
\a) G /C"*- and for all n G {0, 1, 2, . . . , Ai'o — 1}, («|//"|$) = 0. Let P be the projector onto /C and 
let H = PHP be the Hamiltonian within this subspace. Then 
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Proof Define and \^(t)) as 
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Note that by assumption, H^\^) = H''\^) for all k < Nq, and thus the first A'o terms in the two 
above sums are equal. Looking at the difference between these two states, we have 
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where the last step uses the fact that \\H\\ < \\P\\ \\H\\ \\P\\ = \\H\\. Thus for any c > 1, we have 



mt)) - mm < ^-^m' 



c"" k\ 
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<_exp(ct||iJ||). 

We obtain the best bound by choosing c = No/\\Ht\\, which gives 
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as claimed. (If c < 1 then the bound is trivial.) 

Proposition 2. Let Ui, . . . ,Un and Vi, . . . ,Vn be unitary operators. Then for any 



□ 



\i=n i=n 



i=j~l 



(57) 



Proof. The proof is by induction on n. The case n = 1 is obvious. For the induction step, we have 
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where the last step uses the induction hypothesis. 

We are now ready to prove Lemma 2: 
Proof of Lemma 2. For M G N write 
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where in the second line we have used Proposition 2 and in the last step we have used Proposition 1 
and the fact that || |7(t)) || < 1 + S. Now, for some r] > 1, choose 
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rie\\H\\t 
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for < i < T to get 

\\(e-'"' - e-'"')\^)\\ < 2M {6 + r]''''>{l + 6)) 

The choice rj — 2 gives the stated conclusion. □ 

Note that it would be slightly better to take a smaller value of rj. However, this does not 
significantly improve the final result; the above bound is simpler and sufficient for our purposes. 
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